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ABSTRACT 


A  theoretical  and  experimental  investigation  of  mixing  and 
supersonic  combust'on  is  presented.  A  review  of  the  problem  under 
investigation  is  given  in  Section  I.  Then  an  analysis  of  inviscid  flow 
fields  with  a  finite  rate  chemistry  for  a  hydrogen-air  reaction  is 
presented  in  Section  II.  In  Section  III,  an  analysis  of  a  turbulent 
mixing  for  fLows  with  large  density  gradients  with  no  chemical  reaction 
is  presented.  The  analysis  is  compared  with  the  experimental  results. 
In  Section  IV,  the  results  of  the  experiments  in  supersonic  combustion 
are  presented. 
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SECTION  I 


INTRODUCTION 

In  recent  years  the  problem  of  combustion  in  a  stream  which 
moves  at  supersonic  velocity  has  attracted  the  attention  of  several 
groups  of  investigators.  One  of  the  practical  motivations  for  research 
on  this  problem  is  the  hypersonic  ramjet  employing  supersonic  com¬ 
bustion.  In  the  studies  of  the  detailed  processes  invoLved  in  the  releasi 
of  chemical  energy  in  a  supersonic  stream  two  modes  of  combustion 
have  been  considered:  detonations  which  are  controlled  by  a  wave 
mechanism  and  supersonic  diffusion  flames  which  are  controlled  by 
mixing  of  the  fuel  and  oxidizer. 

Consider  briefly  these  two  modes;  in  detonation  the  oxidizer  and 
fuel  are  completely  premixed  upstream  of  the  combustor  in  a  region  of 
the  flow  wherein  the  static  pressure  and  temperature  are  sufficiently 
low  so  that  essentially  no  chemical  reaction  takes  pLace.  When  the 
stream  of  premixed  gases  arrives  at  the  combustor,  it  encounters  wave 
which  produce  an  increase  in  pressure  and  temperature  and  which 
initiate  chemical  reaction.  Initially  the  waves  are  generated  by  the 
geometry  of  the  combustor;  however,  the  wave  pattern  is  directly 
influenced  by  the  heat  release  so  that  under  steady  operating  conditions 
the  geometry,  upstream  flow  conditions  and  heat  release  are  intimately 
reLated.  Indeed,  after  steady  combustion  is  established  the  geometry 
of  the  combustion  could  possibly  be  altered  so  that  a  self-sustaining  wa- 
pattern  is  achieved. 

This  detonation  mode  of  combustion  in  supersonic  flows  has  bee 

studied  by  Gross  and  Chinitz^  Nicholls^  ^ ,  and  more  recently  by 

(1-3-4) 

Rhodes  and  co-workers'  '  ’  .  It  is  reLated  from  a  fundamental  poinl 

of  view  to  detonation  waves  which  have  been  studied  for  roughly  80  year 
(cf.  reference  1-5  for  a  recent,  excellent  review  of  the  research  per¬ 
taining  thereto).  While  this  mode  of  combustion  is  of  considerable 
interest  for  the  investigation  of  the  fundamentals  of  chemicaLly  reacting 

Manuscript  released  23  November  1962  by  the  authors  for  publication 
as  an  ARL  Technical  Documentary  Report. 


I 


flows  in  general  and  of  chemical  kinetics  in  particular,  it  seems  from 
the  available  experimental  results  that  it  will  be  difficult  to  apply  in  an 
operating  engine  system  primarily  because  of  the  aforementioned  inter¬ 
dependence  of  combustor  of  aprnHvnamir  and  heat  release 

effects.  Indeed  operating  limits  associated  with  combustion  instability 
may  be  expected  in  an  engine  system  employing  this  mode. 

The  second  mode  of  combustion,  the  supersonic  diffusion  flame, 
involves  a  fuel  jet  in  a  supersonic  oxidizing  stream,  as  shown  schemati¬ 
cally  inFigurel-I,  The  processes  of  mixing  and  chemical  reaction  are  in 
general  intimately  connected,  the  reactants  at  first  mixing  and  then 
combining  so  as  to  release  chemical  energy.  It  will  be  recognized  that 
this  mode  is  analogous  to  the  fuel  jet  in  low  speed  flow,  apparently  first 
treated  in  1928  by  Burke  and  Schumann^  ^(or  see  reference  1-7).  As 
in  all  viscous  flows  involving  chemical  reaction  it  is  instructive  to  con¬ 
sider  two  characteristic  times,  one  associated  with  diffusion  Tj-j  and  one 
with  the  chemical  reaction,  ~c  .  The  diffusion  time  can  be  associated 
with  the  time  required  by  a  fluid  eLement  with  the  external  stream 
velocity  ue  to  travel  the  distance  in  which  the  veLocity  at  the  axis  is 
cLose  to  that  in  the  external  stream.  The  chemicaL  time  can  be  esti¬ 
mated  from  the  kir.etic  rates  under  the  assumption  of  infinitely  fast 
mixing.  In  the  treatment  of  Low  speed  fuel  jets  it  is  implicitly  assumed 
that  td  »  Tc  so  that  effectively  equilibrium  chemistry  prevails,  i.e,, 
so  that  diffusion  is  rate  controlling.  In  a  practical  supersonic  diffusion 
flame  T j-j  and  Tc  can  be  of  the  same  order  for  some  flow  conditions; 
therefore,  in  the  anaLysis  of  the  combustion  process  the  interplay  of 
chemical  and  fLuid  dynamic  effects  must  be  considered. 

There  are  a  variety  of  fluid  mechanical  and  chemical  problems 
which  must  be  considered  in  achieving  an  understanding  of  the  basic 
features  of  this  mode  of  combustion.  There  has  been  under  way  for 
several  years  research  on  these  problems  at  the  Aerodynamics 
Laboratory  of  the  Polytechnic  Institute  of  BrookLyn  (PIBAL)  and  at  the 
General  Applied  Science  Laboratories  (GASL).  It  is  the  scope  and 
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purpose  of  this  paper  to  review  some  of  the  fundamental  aspects  of  this 
research  effort  including  the  progress  to  date,  and  to  indicate  problem 
areas  reauiring  additional  research. 

For  purposes  of  orientation  it  may  be  instructive  to  consider 
some  of  the  overall  aspects  connected  with  the  hypersonic  ramjet 
employing  supersonic  combustion.  In  the  discussion  of  reference  1-8 
at  the  First  ICAS,  the  first  author  suggested  for  the  first  time  the  use 
of  defLagration  or  diffusion  fLame  in  a  supersonic  combustor.  This 
idea  was  developed  further  in  reference  1-9  where  the  first  author  pre¬ 
sented  a  detailed  discussion  of  a  hypersonic  ramjet  employing  hydrogen 
as  a  fuel  with  this  mode  of  combustion.  It  was  shown  therein  that  such 
a  ramjet  can  produce  a  large  thrust  per  unit  frontal  area  and  can 
operate  at  a  high  specific  impulse  even  at  hypersonic  speeds.  Figure  1-2 
is  taken  from  reference  1-9  and  provides  the  results  significant  for 
present  purposes.  Other  studies  of  various  aspects  of  hypersonic 
ramjets  employing  supersonic  combustion  are  discussed  in  refer¬ 
ences  1-10  and  -11, 

It  is  fundamental  to  the  operation  of  such  a  ramjet  to  decelerate 
the  flow  in  an  inlet  from  the  flight  Mach  number  to  a  supersonic  value 
such  that  at  the  entrance  to  the  mixing  region  the  static  temperatures 
are  in  the  range  of  1000  to  2000°K.  The  hydrogen  fuel  is  injected 
parallel  to,  and  in  the  downstream  direction  with  respect  to,  the  air 
flow.  Mixing  and  heat  reLease  occur  in  the  combustor,  by  controlling 
the  mixing  and  by  changing  the  cross-sectional  area  available  to  the 
flow,  the  pressure  distribution  in  the  mixing  region  can  be  controlled 
so  as  to  avoid  shock  formation.  A  nozzle  connected  to  the  combustor 
expands  the  combustion  products  to  obtain  gross  thrust  from  the  engine. 
In  view  of  the  specific  impulse  which  can  be  realized  with  hydrogen  as 
a  fuel,  research  at  PIDAL  and  GASL,  has  been  confined  to  the  hydrogen- 
air  system;  only  this  system  will  be  discussed  explicitly  herein  although 
many  aspects  of  the  discussion  would  apply  to  other  fuels. 
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FIG.  1-2.  Ramjet,  'nlet  Static  Temperatures ,  Thrust,  and  Specific  Impulse, 
for  Various  Free  Stream  Mach  Numbers  I  Hydrogen  Fuel'. 

Altitude  1  50,000  ft. 
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The  investigation  of  supersonic  diffusion  flames  is  complicated 
by  the  simultaneous  occurrence  of  fluid  mechanical  and  chemical 

s  g »•  o£  £\^wL  oxidiz ° v j  ■*•«•» gnnn >•  n  wi^^i  r ^ p 

heat  release  and,  therefore,  produces  large  changes  of  density  and 
pressure.  In  a  supersonic  combustion  process  the  pressure  variation 
propagates  in  the  flow  through  a  wave  mechanism  in  a  diffused  medium. 
A  complete  analysis  of  such  processes  is  complex  and,  therefore,  many 
simplifications  are  required  to  increase  the  understanding  of  the 
phenomena.  For  a  large  range  of  static  pressure  and  temperatures 
important  for  practical  applications,  the  reaction  process  is  much 
faster  than  the  diffusion  process.  In  this  case  the  effect  of  transport 
properties  during  the  time  in  which  the  chemical  process  takes  place 
can  be  neglected  in  the  first  approximation,  and  the  concentration  of  all 
elements  along  a  given  streamline  can  be  assumed  to  be  constant, 
unaffected  by  the  diffusion  process.  Then  the  chemical  behavior  along 
the  streamlines  can  be  considered  an  inviscid  phenomenon.  Pressure 
and  density  can  change  along  the  streamline;  therefore,  chemical 
reaction  with  variabLe  pressure  must  be  analyzed.  The  results  of  this 
type  of  analysis  will  indicate  for  what  range  of  problems  this  assumption 
can  be  accepted.  If  this  assumption  is  sufficiently  accurate,  the  mixing 
process  can  also  be  analyzed  in  a  simplified  form  by  assuming  that  the 
chemical  reaction  is  infinitely  fast  with  respect  to  the  mixing  process. 
Then  the  flow  everywhere  is  in  chemical  equilibrium;  however,  the 
mixing  may  occur  with  pressure  gradients  in  all  directions.  The 
combustion  process  that  is  of  practical  interest  is  one  which  avoids 
formation  of  shock  waves  in  the  flow;  therefore,  the  analysis  of 
heterogeneous  mixing,  i.e.  ,  of  fuel  and  oxidizer,  with  constant  pres¬ 
sure  is  of  great  importance.  However,  the  mixing  rate  and  heat 
release  must  be  controlled  and  related  to  the  shape  of  the  combustor 
in  order  to  achieve  constant  pressure. 

Following  this  line  of  investigation  this  paper  is  organized  as 
follows:  Section  II  is  devoted  to  the  analysis  of  inviscid  flows  with 
finite  rate  chemistry,  to  a  discussion  of  the  mechanisms  in  the 


hydrogen-air  reaction,  and  to  the  description  of  the  relation  between 
combustor  shape  and  mixing  process  required  to  avoid  shock  formation. 
Section  III  is  concerned  with  turbuLent  mixing  in  flows  with  large  density 
gradients  due  to  large  differences  in  temperature  and  in  concentration 
of  Light  gases  but  with  either  no  chemical  reaction  or  with  equilibrium 
chemical  behavior  Finally,  the  results  of  experiments,  which  employ 
the  hypersonic  faciLity  of  PIBAL  to  investigate  the  combustion  of 
hydrogen  in  a  supersonic  stream,  are  discussed  in  Section  IV. 

The  authors  gratefully  acknowledge  the  support  and  cooperation 
of  the  Air  Force  organizations  and  the  efforts  of  Mr,  Samuel  Lederman, 
Mr.  Louis  ALpinieri,  Mr.  GdaLia  Kleinstein,  and  Mr,  Herbert  Fox,  and 
of  Mr.  Harold  Pergament,  In  addition,  the  authors  acknowledge  with 
thanks  helpfuL  discussions  with  Professors  ManLio  Abele  and  Lu  Ting. 


SECTION  II 

INVISCID  FLOWS  WITH  CHEMICAL  REACTION 

As  part  of  the  investigation  of  supersonic  diffusion  flames  with 
the  hydrogen-air  system,  there  is  considered  herein  the  mechanisms 
and  kinetics  involved  when  mixing  and  diffusive  effects  are  absent.  These 
considerations  correspond  to  the  assumption,  either  of  premixed  fuel  and 
oxidizer,  or  to  the  other  limit  wherein  chemical  reaction  is  so  fast  with 
respect  to  the  diffusion  process,  so  that  the  diffusive  effects  taking  place 
during  the  combustion  process  can  be  neglected. 

The  problem  of  inviscid  gas  flows  with  non-equilibrium  chemical 
behavior  has  been  studied  intensively  in  the  past  several  years.  Primary 
attention  has  been  devoted  to  non-equilibrium  air  flows  around  bodies  in 
hypersonic  flight  and  through  nozzles  of  hypersonic  test  facilities.  Recently, 
nozzle  flows  involving  combustion  products  of  interest  in  high-performance 
propulsion  devices  have  been  considered  because  of  the  practically  important 
differences  in  gross-thrust  associated  with  the  limiting  cases  of  frozen  and 
equilibrium  flow.  Olson ^  ^  provides  a  recent  review  of  these  latter  fLows. 

In  all  of  the  analyses  of  non-equilibrium  behavior  for  propulsion 
the  flow  fieLds  and  the  thermodynamic  behavior  of  the  gases  have  been 
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idealized.  The  flow  is  treated  either  as  one-dimensional  with  the  area- 
distribution  of  the  nozzLe  specified  or  is  considered  aLong  streamlines 
with  the  pressure  distribution  specified.  As  wiLl  be  discussed  in  detail 
beLow,  the  latter  idealization  appears  to  have  several  advantages.  With 
respect  to  the  thermodynamic  behavior  of  the  combustion  products, 
there  will  be  assumed  thermodynamic  equilibrium,  i.e.,  full  rotationaL 
and  vibrational  equilibration  of  all  molecules  at  a  local  t ransLational 
temperature  which  is  the  same  for  all  species.  In  addition,  the  forward 
and  reverse  rate  constants  are  assumed  to  be  related  by  the  appropriate 
equilibrium  constant.  It  is  generally  considered  that  these  assumptions 
are  satisfactory  in  view  of  the  inaccuracy  with  which  the  mechanisms 
and  rates  associated  with  practical  fuel-oxidizer  systems  are  known. 

1.  REVIEW  OF  RECENT  RESEARCH 

Consider  now  the  available  results  of  theoretical  calculations 
applicable  to  the  hydrogen-air  system.  The  fundamental  paper  is  that 
of  Duff^  ^  who  carried  out  calculations  of  the  reaction  profile  behind 
a  detonation  wave  satisfying  the  Chapman- Jouget  condition  in  a  hydrogen- 
oxygen  mixture.  Duff  attributes  the  individual  reaction  steps  and  the 
reLevant  rate  data  to  N.  Davidson.  This  research  was  continued  at 
Los  Alamos  by  Schott,  et  al.  experimental  research  confirmed 

in  generaL  the  overalL  accuracy  of  the  chemistry.  In  reference  II-5, 
Libby,  Pergament,  and  Bloom  presented  the  results  of  reaction  histories 
for  hydrogen-air  mixtures  at  constant  pressure  using  the  chemistry  of 
Duff  and  nitrogen  treated  as  an  inert  diluent.  Recently,  Westenberg 
and  Favin^  ^  carried  out  one-dimensional  nozzle  calculations  for  the 
hydrogen-air  system  and  for  a  second  mixture  corresponding  to  rocket 
exhaust  gases.  A  range  of  pressures  in  the  settling  chamber,  a  single 
stagnation  temperature  of  3000°K  and  two  equivalence  ratios  were 


Some  of  the  results  of  these  and  similar  calculations  will  be  discussed 
in  detail  below. 
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assumed  There  was  obtained  the  expected  effect  of  pressure,  namely, 
that  the  flow  tends  to  remain  near  equilibrium  when  the  pressure  is 
high  and  to  depart  from  equilibrium,  indeed  to  freeze  effectively  at  some 
point  downstream  of  the  throat,  when  the  pressure  is  low.  Also  of 
interest  "'is  a  result,  found  independently  by  Y/c stenbe r g  and  IT avm  and 
by  the  present  authors  as  discussed  below,  concerning  the  behavior  of 
the  intermediates  atomic  hydrogen  and  the  hydroxyl  radical  under 
essentially  frozen  conditions.  It  is  found  that  in  an  expanding  flow  at 
sufficiently  low  pressures  so  that  the  three-body  reactions  are  frozen, 
the  intermediates  do  not  necessarily  have  concentrations  between  their 
frozen  and  equilibrium  values. 


Finally,  Momtchiloff,  et  al.  ®  ^  carried  out  calculations  for 
constant  area  combustion  of  hydrogen  and  air.  A  range  of  initial  con¬ 
ditions  of  pressure,  temperature,  velocity  and  equivalence  ratio  were 
covered  although  to  reduce  computing  time,  the  reaction  profiles  were 
not  carried  to  equilibrium.  In  this  work  several  nitrogen-oxygen 
reactions  leading  to  the  formation  of  nitric  oxide  (NO)  were  included 
in  the  reaction  steps  but  were  found  to  be  unimportant  below  2200°K, 
justifying  in  part  the  assumption  employed  herein  of  nitrogen  as  a 
diluent. 


As  discussed  in  more  detaiL  below,  the  reaction  steps  and  the 

rate  constants  for  the  above  references  are  essentially  the  same  as 

those  given  by  reference  II- 2.  The  rate  constants  differ  somewhat  but 

>!< 

probably  only  within  the  accuracy  to  which  such  data  are  known. 


After  completion  of  this  manuscript  the  attention  of  the  authors  was 
called  to  the  recent  paper  by  Fowler.  R.  G.  ,  "A  Theoretical  Study  of 
the  Hydrogen -Air  Reaction  for  Application  to  the  Field  of  Supersonic 
Combustion."  Proceedings  of  the  1962  Heat  Transfer  and  Fluid 
Mechanics  Institute,  Stanford  University  Press.  There  is  presented 
constant  density  calculations  employing  chemical  mechanisms  and 
rates  similar  to  those  described  therein. 


9 


2.  ANALYSIS  ALONG  STREAMLINES  -  GENERAL  CONSIDERATIONS 


To  indicate  in  more  detail,  the  characteristics  of  the  hydrogen- 
air  reaction  in  a  flowing  system  and  the  methods  of  anaLysis  for  inviscid 
flows,  the  pertinent  gas  dynamical  equations,  tie  reaction  steps,  the 
rate  constants,  and  some  numerical  results  will  be  discussed.  Con¬ 
sider  the  inviscid,  steady  fLow  along  streamlines  of  a  chemical  reacting 
gas  involving  N  species  and  L  elements.  The  conservation  of 
momentum,  energy,  species  and  elements  can  be  stated  as  follows: 


Momentum 


Energy 


pV  f  f 

ds  ds 


N 


d H 
ds 


=  0,  H  =  {V2/  2)  +  ;  Y.h. 


Species 


Elements 


dY. 

pV  - -  =  w.  ,  i  =  1,  2,  •  •  •  N-L 

ds  i 


dY. 


ds 


=  0  ,  j  =  N-L  +  1,  ...  N 


N 


Y.  s  )  u-.W.Y./W. 

J  L  MiJ  J  x'  i 

i=l 


(H-2.1) 

(II- 2.2) 

(II- 2.3) 


(II- 2.4) 


These  equations  are  supplemented  by  algebraic  equations:  An  equation 
of  state 


N 

p  =  0  R  T  )  Y./W.  =  pR  T/W  ;  (II- 2.5) 

o  1  1  o 

i=l 
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the  creation  terms  as  functions  of  composition  and  state;  and  the  species 
enthalpy-temperature  relations  in=h.{T),’1' 

The  form  of  the  creation  terms  is 

K  N 

W  =  W  )  (v'  -  v/.)k  onj  TT  (Y  /W  (II- 2.6) 

1  1  f-»  lJ  J  .  7  11  J 

j=l 

where 

m.  N 

G.  h  1  -  0  J  K;  TT  (Y./W  )  lJ  1J  (II- 2.7) 

J  c,j  ,=1  i  i 

mj=I  (viTV 

i=i 

and  where  k.  is  the  specific  raue  constant,  K  .  is  the  equilibrium 
J  C  j  J 

constant  based  on  molar  concentrations,  and  vf.  and  sj.'.  are  the 

th  G  1J 

stoichiometric  coefficients  for  the  j  reaction 

N  k  N 

l  vA  M. - ^  vfj  M.  ;  j=l,  •••  K  (n-2.8) 

i=I  i=l 

Basically,  these  equations  have  been  employed  to  date  in  all 
inviscid  calculations  with  non-equilibrium  chemical  behavior.  In 
nozzle  flows,  the  system  of  equations  is  usually  completed  by  the 
assumption  of  one-dimensional  flow  so  that  pVA  =m  with  the  area 
distribution  prescribed. 

”*  *“  *  “*  *  '  1  ' — *  ' — ' — '  * 1 

In  the  calculations  reported  here  linear  relations  between  h^  and  T  have 
been  assumed.  This  corresponds  to  taking  an  average  Cp  {  for  the 
temperature  range  of  interest  cind  leads  through  Eq.  (II- 2/2)  to  a  con¬ 
venient,  explicit  relation  for  temperature  T  =  T (Y.  ,  V3) . 
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This  approach  does  not  permit  generalization  to  more  complex 
flow  field  calculations,  for  example,  performed  in  connection  with  the 
method  of  characteristics.  In  such  analyses  the  physical  properties 
such  as  Pressure  and  density  are  determined  aLong  any  given  Stream- 
Line,  and  not  streamtube  area.  For  these  reasons  herein  the  pressure 
distribution  is  specified  along  the  streamline  under  consideration,  as 
suggested  by  Bloom  and  Steiger ^  ^  in  connection  with  non-equilibrium 
dissociation  and  recombination  in  external,  hypersonic  flows.  This 
also  completes  specification  of  the  problem,  i.e.  ,  given  initial  condi¬ 
tions  at  a  generic  point  (s  =  0)  the  distribution  of  fLow  composition  and 
state  variables  with  respect  to  s  can  be  found  by  integration. 

It  is  perhaps  worth  while  at  this  point  to  make  some  remarks 
about  several  numerical  problems  which  arise  in  the  treatment  of  the 
above  equations.  The  number  of  differential  equations  which  must  be 
integrated  with  respect  to  the  space  variable  s  can  be  reduced  by 
employing  the  element  conservation  equations  [Eqs.  (II- 2.4) ].  These 
can  be  integrated  immediately  and  provide  L  algebraic  equations  among 
the  N  species,  so  that  only  N-L  equations  involving  creation  terms 
need  be  considered.  However,  in  many  calculations  the  mass  fractions 
of  the  intermediates  can  be  many  orders  of  magnitudes  Less  than  the 
mass  fractions  of  the  major  constituents.  As  a  result  the  solution  of 
the  algebraic  equations  can  invoLve  a  loss  of  accuracy  due  to  the 
practical  problem  of  small  differences  in  relatively  Large  numbers. 
Similar  difficulties  can  be  encountered  when  the  flow  is  near-equilibrium; 
the  creation  term  in  Eqs.  (II-2.3)  represents  in  generaL  the  sum  of 
differences  between  the  forward  and  backward  rates  associated  with  a 
series  of  reaction  steps.  At  equilibrium  in  allowing  system  with 
pressure  and  temperature  changes,  the  net  result  of  the  sum  is  non¬ 
zero  but  can  invoLve  smaLl  differences  and  thus  a  loss  in  accuracy. 
Another  numerical  probLcm  is  rcLatcd  to  starting  the  integration  in 
some  cases;  in  nozzle  calculations  the  combustion  products  in  the  settling 
chamber  or  combustor  are  usuaLly  assumed  to  be  in  equilibrium.;  but  at 
equilibrium  the  creation  terms  must  be  considered  indeterminant  in  that 
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a  Large  reaction  rate  multiplies  a  quantity  which  is  the  equilibrium 
condition,  i.  e.  ,  a  quantity  which  is  essentially  zero.  The  product 
yieLds  a  value  of  \i\  maintaining  the  species  concentration  at  its  equi¬ 
librium  value  at  each  point  in  the  flow.  Thus  integration  from  an 
equilibrium  state  invoLves  a  singular  perturbation  such  as  has  been 
discussed  by  BLoom  and  Ting^  ^  .  HalL  and  Russo®  and 
Vincenti^  However,  it  has  been  found  possible  in  actual  calcu¬ 

lations  as  per  the  present  report  and  in  references  TI-6  and  -12  to  start 
the  integration  from  a  nominal  equilibrium  state  with  an  initiaL  step- 
size  which  is  smaller  than  that  usually  required  by  integration  accuracy. 
After  a  few  integration  steps,  the  effect  of  the  numerically  uncontrolled 
start  is  found  to  disappear  and  the  solution  of  the  differential,  equations 
becomes  a  welL-behaved  approximation  to  equilibrium  flow,  which  is,  of 
course,  described  by  a  solution  to  a  set  of  algebraic  equations.  Finally, 
it  is  remarked  that  computing  time  with  standard  accuracy  requirements 
for  the  selection  of  step-size  can  become  excessive  even  for  high  speed 
computing  equipment  when  the  flow  is  near-equilibrium.  This  is  due  to 

3jc 

a  tendency  for  the  solution  to  oscillate  about  the  equilibrium  solution, 

3.  THE  REACTION  STEPS  AND  REACTION  RATES 

The  reaction  steps  and  rate  constants  used  here  are  derived 
from  the  work  of  Duff  and  co-workers^  ^  t°  4)^  same  chemical 

reactions  suggested  by  Duff  have  been  used  independently  by  the  several 
investigators  cited  above;  thus  it  currently  appears  to  represent  the 
best  available  set  of  reactions  for  the  hydrogen-oxygen  system.  Here 
and  in  reference  H-6  nitrogen  is  treated  as  an  inert  diluent;  this 
assumption  as  shown  in  reference  IL-7  is  considered  valid  for  temper¬ 
atures  beLow  roughLy  2500°K.  It  is  worth  mentioning  at  this  point  that 
the  temperature  range  of  interest  in  supersonic  combustion  appears  to 


A  detailed  discussion  of  numerical  difficulties  associated  with 
caLculatimis^qf  non-equilibrium  nozzle  flows  is  given  by  Emanuel  and 
Vincenti'  ',  The  current  state  of  development  for  air  calculations 
is  indicated  in  reference  11-13, 
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be  below  roughly  3000°K;  this  is  within  the  range  covered  by  the  shock 
tube  experiments  of  Schott  et  al.  ®  ^ ^  and  only  somewhat  above 
the  temperature  obtained  in  detonation  experiments,  e.g.  ,  Nicholls®  ^ 
and  in  laminar  flame  experiments,  e.g.  ,  Fine®  These  experi¬ 

ments  do  not  verify  and  establish  the  individual  reaction  steps  and 
reaction  rates  but  have  verified  partially  the  overall  consistency  of 
the  ensemble  required  for  practical  analysis. 


The  reaction  rates  considered  here  are  as  follows: 

kf 


1 

H  +  Oa 

- 

OH 

+  O 

.  -8810/T 

3(1 01 4 )  e 

2 

o  +  h3 

-« 

OH 

+  H 

, „  - 4030/T 

3(1 01 4 )  e 

3 

OH  h  H2 

- 

h3° 

+  H 

3(10>.).-3020/T 

4 

20  H 

— ♦ 

Ha° 

+  O 

3(10-)e-3020/T 

RR  No.  1  No.  2 

No.  3 

5 

2H  +  M 

— « 

Hs 

+  M 

10ls  5(1015) 

K101  ®) 

6 

H  +  OH  + 

M  - 

h2o 

+  M 

1016  1017 

3(1017) 

7 

H  +  O  + 

M  - 

OH 

+  M 

1016  10lc 

3(I016) 

8 

20  +  M 

- 

o3 

+  M 

3(1C14) 

The  temperatures 

are 

in  °K;  the  forward  rate  constants. 

i.  e.  ,  fo 

(II- 3.1) 


-l  -i 

sec  for 


reaction  proceeding  to  the  right  are  given  in  (mole/cc) 

_  2  _  l 

second  order  reactions  (1-4)  and  in  (mole/cc)  sec  for  third  order 
reactions  (5-8);  the  symbol  M  is  used  to  denote  any  molecule  or  atom 
serving  as  a  third  body.  As  will  be  discussed  in  detail  below,  the 
rate  constants  for  reaction  steps  5  and  6  play  an  important  role  at  low 
pressures  in  determining  the  chemical  history  of  the  flow;  therefore, 
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several  values  thereof  have  been  selected.  Those  denoted  RR  No.  1 
and  RR  No.  2  represent  the  range  of  probable  values  estimated  by 
Schott^  ^  while  RR  No.  3  are  values  considered  possible  on  the  basis 
of  more  recent  combustion  results.  The  rate  constants  for  the  reverse 
reactions  are  obtained  from  the  equilibrium  constants  approximated  as 
functions  of  the  temperature  in  the  form  Be^/^. 

To  indicate  the  range  of  values  of  the  reaction  rate  constants 
employed  in  related  theoretical  studies  there  are  compared  in  Table  II- 1 
the  rate  constants  for  reaction  5  employed  therein;  this  reaction  is 
selected  as  being  an  important  but  typical  one  in  the  array  of  reaction 
steps.  It  wilL  be  noted  that  in  references  II-3,  -6,  and  -7  a  different 
reLation  between  the  temperature  and  ks  has  been  assumed.  However 
in  the  temperature  range  of  1000  to  3000°K  the  numerical  values  given 
by  the  different  expressions  are  within  an  order  of  magnitude;  this  is 
indicative  of  the  accuracy  with  which  these  physical-chemical  data  are 
known. 


TABLE  II-l 

Reaction  Rates  for  Hydrogen  Recombination 


2H  +  M  -  H2  +  M 


Present  Report  RR  No.  1 
RR  No.  2 
RR  No.  3 


Reference  H-3 
Reference  II-6 
Reference  H-7 


ks(moles/cc)  sec 

101  B 
5  ( 1 0 1  s ) 

1.2(1016) 

10ls  to  5 (10x  6) 

2(I018)T_1 

5(1016)T‘2 


The  reaction  steps  described  by  Eqs.  (II -3.1)  can  be  separated 
into  two  groups;  the  second-order  or  "shuffling"  reactions  (1-4)  which 
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proceeding  to  the  right  deplete  the  population  of  molecular  oxygen  and 
hydrogen  and  lead  to  the  formation  of  water  and  of  the  intermediates, 
atomic  oxygen  and  hydrogen  and  the  hydorxyl  radical;  and  the  third 
order  reactions  (5-8)  which  proceeding  to  the  right  deplete  the  inter¬ 
mediates  and  lead  to  the  formation  of  water  and  molecular  species. 

This  complex  array  of  forward  and  reverse  reactions  makes  difficult, 
if  not  impossible.,  a  priori  estimates  of  chemical  behavior  in  general 
and  of  chemical  times  in  particular. 

4,  NUMERICAL  RESULTS  FOR  CONSTANT  PRESSURE 

It  is  instructive  for  an  understanding  of  the  kinetics  involved 
in  the  supersonic  combustion  of  hydrogen  to  consider  the  reaction 
history  of  the  following  flow  model:  a  mixture  of  hydrogen  and  air  at 
time  zero  at  an  initial  temperature  of  1000°K  or  higher.  The  pressure 
is  assumed  constant  with  time. 

In  the  calculations  presented  here,  some  of  these  results  are 
taken  from  reference  II-5;  the  additional  calculations  have  been  per¬ 
formed  by  H.  Pergament  and  will  be  reported  "in  extenso"  elsewhere, 
the  initial  state  of  the  gas  is  assumed  to  correspond  to  a  partial 
equilibrium  between  the  atomic  and  molecular  species  of  oxygen  and 
between  the  atomic  and  molecular  species  of  hydrogen.  The  results  are 
not  significantly  aLtered  if  the  initial  concentrations  of  the  atomic  species 
are  changed  somewhat;  relatively  large  concentrations  of  atomic  hydrogen 
would  of  course  change  the  early  phases  of  the  reaction  histories. 

Consider  the  temperature  histories  for  some  typical  cases 
based  on  this  model  as  shown  in  iigures )'  -  1j.-c  and  as  computed  by  a 
standard  integration  program  on  an  IBM  7090  computer.  Presented  are 
the  results  for  three  pressures,  for  several  equivalence  ratios’’  and  for 
the  three  rates  of  reactions  5  and  6  as  discussed  in  Section  II-3  above. 


The  equivalence  ratio,  denoted  as  -n  ,  is  defined  by  the  chemical  equa¬ 
tion  for  the  reactants;  namely..  +  2-^H3  +  3-76N3  -•  products. 
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Temperature-time  Histories  of  Hydrogen-Air  Mixtures  at  Constant  Pressur 


=  1  ATM 


TIME 


TIME  -  SECONDS 


FIG.  II-lc.  Temperature-Time  Histories  of  Hydrogen-Air  Mixtures 
at  Constant  Pressure. 
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These  temperature  histories  are  characterized  by  two  time  periods; 

there  is  an  initial,  time  interval,  usually  termed  the  induction  time  and 

denoted  t.  ,  during  which  the  temperature  remains  sensibly  constant; 

and  a  later  phase,  during  which  the  temperature  increases  from  its 

initial  value  T.  to  the  final  temperature  T^.  the  latter  corresponds  in 

this  model  to  the  adiabatic  flame  temperature  (reference  11-17).  The 

time  associated  with  the  second  phase  is  termed  here  the  energy  release 

time  and  is  denoted  t  :  thus  r  =  t  +  t  .  Under  conditions  of  high 
r  c  1  r 

initial  temperature  and  low  pressure  the  second  time  interval  dominates 
while  under  conditions  of  relatively  low  initial  temperature  {I0QQ°K) 
and  high  pressure  (4  atmospheres)  the  two  time  intervals  axe  of  roughLy 
equal  duration.  This  last  result  appears  not  to  have  been  observed 
previously. 

The  changes  in  composition  during  these  periods  are  shown  in 
representative  cases  in  Ilgurcsj  J- 2a-f.  Of  greatest  thermodynamic 
significance  because  of  its  high  species  enthalpy  per  unit  mass  is  the 
behavior  of  the  atomic  hydrogen'  during  the  induction  phase  the  con¬ 
centration  of  H  and  the  other  intermediates  O  and  OH  increases  rapidly 
and  reaches  a  maximum  corresponding  to  mass  fractions  of  the  order 
of  several  percent  or  less.  At  the  same  time  the  molecular  species, 

Os  and  Hs  are  depleted  and  water  is  formed.  However,  the  chemical 
energy  reLeased  by  the  formation  of  water  is  effectively  absorbed  by 
the  atomic  hydrogen  so  that  no  significant  change  in  temperature  occurs. 
During  the  second  time  period  the  intermediates  decay  to  their  appro¬ 
priate  equilibrium  vaLues  releasing  chemical  energy  and  leading  to  the 
temperature  rise  seen  in>  inures  J.i  -  la -c.  It  should  be  noted  that  at  the 
end  of  the  induction  period  the  major  constituents  are  close  to  their 
equilibrium  values;  during  the  second  period  the  decay  of  the  inter¬ 
mediates  leads  to  only  secondary  adjustments  m  the  concentrations  of 
the  major  constituents. 

This  behavior  is  explicable  in  terms  of  the  reaction  steps  given 
by  Eqs  (JI-3.1)  During  the  induction  phase  the  two-body  or  "shuffling" 


20 


21 


FIG.  II- 2a.  Composition- 


FIG 
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FIG.  II-?.d.  Composltion-  TIrm  Hl-'tories  of  Hydrogen-Air  Mixtures 


FIG.  n-2f.  C  ompo8ition-Time  Histories  of  Hyd  rogen-Ai  r 
Mixtures  at  Constant  Pressure,  HO. 
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reactions  are  relatively  fast  and  lead  to  an  overproduction  of  the  inter¬ 
mediates;  during  the  second  period,  the  third-order  reactions  lead  to 
recombination  and  depletion  of  the  intermediates  with  the  accompanying 
energy  release.  This  characteristics  of  the  hydrogen-oxygen  reaction 
was  apparently  published  first  by  Duff®-^.  The  results  presented 
here  confirm  this  behavior  over  a  variety  of  initial  conditions. 

It  is  of  interest  to  note  that  at  Low  pressures  the  overalL  reaction 
time  is  sensitive  to  the  recombination  reaction  rate,  e.g.,  for  p=0.2 
atmospheres  the  overalL  reaction  time  is  15  times  Less  for  RR  No.  3 
than  for  RR  No.  1.  On  the  contrary,  for  p=  4.7  atmospheres  these  times 
differ  by  only  a  factor  of  three.  It  must  also  be  noted  that  for  static 
pressure  of  the  order  of  one  atmosphere  or  higher,  the  energy  reLease 
time  is  short;  therefore  the  difference  between  the  resuLts  for  different 
reaction  rates  is  not  important  for  practical  applications.  However, 
this  is  not  the  case  at  Low  pressures. 

5.  THE  INDUCTION  PERIOD 

In  the  past  considerable  attention  has  been  devoted  to  the  time 
associated  with  the  induction  period;  for  example,  Schott  and  Kinsey^  ^ 
and  NichoLLs^  present  the  results  of  experiments,  employing 
respectively  techniques  associated  with  the  shock  tube  and  standing 
detonation  waves,  to  show  the  dependence  of  T.  on  pressure,  temper¬ 
ature  and  composition;  in  both  reference  theoretical  considerations  based 
on  the  reaction  steps  were  employed  to  facilitate  correlation  of  the  data. 

A  convenient  formula  which  can  be  used  for  estimating  T,.  within 
a  factor  of  two  or  three  for  mixtures  with  equivalence  ratios  on  the  order 
of  unity  can  be  obtained  by  simplifying  the  functional  form  suggested  by 
NichoLLs  and  by  employing  coefficients  selected  from  the  present  numerical 
resuLts.  With  t.  defined  as  the  time  at  which  the  computed  temperature 
starts  to  rise  significantly,  there  is  obtained 
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(II  -5.1) 


4  1  04/t  . 

PTX  ~  45(10"  ) e  1 

where  is  in  |j  sec,  p  in  atmospheres  and  in  °K .  In  Figure  II - 3  the 
numerical  results  are  compared  with  the  predictions  of  Eq.  (II-5.1);  it 
will  be  seen  that  within  the  cited  factor  of  two  or  three  the  correlation 
is  satisfactory. 

From  Figures  II-la-c  and  from  11-3  it  will  be  noted  that  the 
induction  time  is  not  greatly  altered  by  changes  in  the  equivalence  ratio 
nor  in  the  reaction  rates  for  recombination.  By  assuming  the  third 
order  reactions  to  have  a  negligible  role  and  the  temperature  to  be 
constant  during  the  induction  period,  it  is  possibLe  to  make  estimates 
of  the  composition  at  the  end  of  the  induction  period.  Such  estimates 
were  made  by  Schott^  ^  and  were  shown  to  be  in  good  agreement  with 
more  complete  calculations  and  with  shock  tube  experiments.  The 
calculation  assumes  that  at  the  end  of  the  induction  phase  there  exists 
a  partial  equilibrium  state  corresponding  to  equilibrium  of  any  three 
of  the  two-body  reactions  (1-4).  In  addition,  since  these  dominant 
reactions  involve  no  change  in  moles  of  gas,  and  since  the  temperature 
is  approximately  constant,  it  can  be  assumed  that  the  molecular  weight 
of  the  mixture  remains  constant  during  the  induction  period.  The 
resulting  three  equilibrium  conditions,  the  statement  concerning 
moLecular  weight,  and  element  conservation  determine  the  mass 
fractions  of  aLl  species.  The  parameters  involved  in  such  a  calculation 
are  r|  and  T^  ;  the  pressure  does  not  enter. 

In  Figure  H-4  the  mass  fractions  of  atomic  hydrogen  at  the  end 
of  the  induction  period  as  computed  by  this  scheme  are  shown  for  a 
range  of  initial  temperatures.  The  calculations  were  carried  out  on  a 
BendixG-15  computer  and  involved  solution  by  a  combination  of  iteration 
and  Newton- Raphson  methods  of  a  set  of  algebraic,  nonlinear  equations. 
Mr.  H.  Fox  carried  out  these  calculations.  Also  shown  are  the 
maximum  values  of  Yc  obtained  from  the  computer  program.  The 
satisfactory  agreement  will  be  noted;  similar  agreement  is  found  for 
the  remaining  constituent. 
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FIG.  H-3.  Cor  relation  of  Compiled  Induction  Times 


t>.  APPLICATION  OF  PARTIAL  EQUILIBRIUM  FOR  ESTIMATES  OF 
THE  ENERGY  RELEASE  PERIOD 

As  a  result  of  the  estimates  of  the  induction  time  and  of  the 
scheme  described  above  for  approximating  the  composition  at  the  end 
of  the  reaction  period,  it  is  possible  to  consider  a  further  approxi¬ 
mation  applicable  to  the  energy  release  period.  The  basis  for  the 
approximation  resides  in  the  coexistence  during  this  period  of  second 
and  third  order  reactions,  i.e.  ,  those  corresponding  to  reactions  1-4 
and  to  5-8  ,  respectively;  if  it  is  assumed  that  the  former  are  suffi¬ 
ciently  fast  relative  to  the  latter  so  as  to  maintain  during  the  entire 
second  period  a  partial  equilibrium  state,  an  approximate  description 
of  the  reaction  history  can  be  obtained  by  a  single  quadrature.  It  is 
noted  that  the  concept  of  partial  equilibrium  can  be  applied  onLy  when  the 
pressure  is  sufficiently  high  so  that  the  assumed  equilibration  of  the 
two-body  reactions  can  occur.  Analyses  of  this  type  have  been  carried 
out  by  Westenberg  and  Favin^  ^  and  by  Zinman  and  Romano®''^. 

Care  must  be  exercised  in  connection  with  the  treatment  of  partial 
equilibrium  in  flowing  systems  in  that  there  are  introduced  mdetermi- 
nacies  which  correspond  to  an  infinite  reaction  rate  multiplying  a  zero 
corresponding  to  the  equilibrium  condition  for  the  fast  reaction  step. 

The  second  author  discussed  this  difficulty  in  general  in  reference  11-19 
and  showed  how  it  may  be  overcome. 

Consider  the  analysis  of  the  second  period  according  to  a  partial 
equilibrium  concept.  Take  for  the  three  independent,  two-body  reactions 
those  denoted  1,  Z,  and  3  in  Eqs.  (II -3.1)  and  take  reactions  5  and  6  to  be 
the  most  important  three-body  reactions ,v  Now  the  assumption  that  the 
former  are  relatively  fast  implies  from  Eq.  (II -2.7) 

On  *1  -  (W„W1  /W7Wb)K- '-.(Y^AeYjl-O  (II- 6,1 

C  )  - 


The  analysis  is  not  complicated  in  principle  by  the  inclusic-i  of  more 
three -bodv  reactions  although  numerically  the  complexity  '  *:  incr6dcpd 
thereby. 
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Gg  =  1  -  (WBW2/W7WG)K^a  (Y,Y6/YbYs)  -  0 
G3  -  1  -  (W7Wa/W3Ws)K;;3  (Y3Y6/Y7Ys)  =  0 


(II- 6.2) 


(H-6.3) 


The  element  mass  fractions  for  the  oxygen-hydrogen  system  with 
nitrogen  treated  as  a  diluent  are  from  Eqs.  (II-2.4) 


Y  =  Y  +  (W  /2W  ) Y  +  Y  +  (W  / 2W  )Y  =  constant  (II- 6.4) 

Y3  =  Ya  +  (W3/W3)Y3  +  Ye  +  (W2/2W7)Y,  =  constant  (n-6.5) 

Y4  =  constant  (II- 6. 6) 

Thus  Eqs.  (II- 6.1)  -  (II- 6.6)  represent  six  equations  in  seven  mass 
fractions  and  in  the  temperature  and  must  be  supplemented  by  a  single 
reaction  equation  and  by  an  energy  equation.  Since  the  atomic  hydrogen 
is  thermodynamically  the  most  important  intermediate,  consider  the 
species  conservation  thereof;  prior  to  making  the  assumption  of  partial 
equilibrium,  i.  e.  ,  thatGk~0,  k=l,  2,  3,  Eqs,  (II—  2.3) ,  (H-2.6)  and 
(II- 2.7)  for  i=6,  give 


dY 

6 

dt 


w 

6 


w 


(W6)  = -  [-k1(Y6/W8)(Yl/Wi)p3Gl+ks(YB/WB)(Y9/Wa)p3Gg  + 


+  k3(Y7/W7)(Ys/Wa)p3G3  -  2k  (Y  /W  )8  (p3/W)G 


-  ke  (Ye  /W6 )  (Y7  /W7 )  (p3  /W)G6  ]  =  F' 


(H-6.7) 


N 


v. . 


Now  with 


partial  equilibrium  assumed,  the  products  k.p  ^  1  T  (Y./W.)  ^G. 

J  i=l  J 

must  be  interpreted  as  an  infinite  rate  k.  multiplied  by  zero,  namely  G-, 

J  J 

such  that  partial  equilibrium  prevails.  Accordingly,  let  cn  denote  these 

terms  and  consider  a.  ,  i  =  1.  2.  3  as  unknowns  but  a*  and  oa  as 
J  "  “ 
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determinant  values.  If  —  0,  k=l,  2,  3  then  Eqs.  (II- 6.1) - (II- 6.3) 
upon  total  differentiation  imply 


,  dT,  8Gk 


>  —  w  k=l,  2,  3 

L,  av  i 
i=l  “*i 


(II- 6.8) 


dT  ,  d  ,  t,  , 

=  P  dT  ( dT  ^nKc,k> 


At  this  point  it  is  necessary  to  consider  the  energy  equation;  for  con¬ 
stant  pressure  flow  the  energy  equation  is  simply  h  =  constant  so  that 


»  w  -  - 1  I  hh’/q 

i=l 


(II- 6.9) 


In  Eqs.  (II-6.8)  and  (H-6.9),  can  be  expressed  in  terms  of  the  cn's 
according  to 


w.  =  W.  Y  (v.".-v/.)a.  ,  j  =  I, 2, 3, 5, 6 
i  i  L  ij  ij  j  J 


(n-6.10) 


so  that  Eqs,  (II-6.8)  and  (II- 6.9)  yield 


—  -  8G,  h.  , 

U<iW  -r  -k  itlKc,k^ViVViVajWi  =  »•  J  =  ‘■2-3-5-6 


L  L  'ay.  c  dT  c»k' '  Ij  lj'  j  1 
i-1  j  1  P 


k  =  1,2,3 


(II-6.11) 


Eqs.  (II- 6.1 1)  permit  ou  ,  j  =  1 , 2, 3  to  be  expressed  in  terms  of  a  $  and 
Os  so  that  the  right-hand  side  of  Eq.  (11-6.7)  involves  only  determinant 
values.  Moreover,  Eq.  (II- 6.7)  can  be  integrated  to  yield 


=  Ti  +  I 


(n-6.12) 
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where  the  constant  of  integration  is  selected  so  that  when  t  -  l\  ,  Ye 
is  equaL  to  that  given  by  the  anaLysis  of  the  induction  period,  namely 

Y®  f  i 

The  calculations  involved  in  applying  this  analysis  can  be 
carried  out  on  a  small  scale  computer;  e.  g.  ,  the  results  reported  here 
were  obtained  on  a  Bendix  G-15  computer.  The  procedure  is  as 
follows:  The  vaLue  of  Ye  varies  from  the  maximum  computed  according 
to  partial  equilibrium  applied  to  the  induction  period  to  that  prevailing 
at  equilibrium  at  the  final  temperature  T£  .  This  latter  value  can  be 
estimated  by  an  equilibrium  calculation  or  from  available  tables  of  the 
combustion  products  of  hydrogen  and  air.  This  range  of  values  is 
divided  into  a  number  of  intervals  appropriate  to  the  integration  of 
Eq.  (II- 6.1 2).  With  a  generic  value  of  Ys  selected  Eqs.  (II- 6.1) - (II- 6.6) 
and  the  energy  equation  h=  constant  are  employed  to  determine,  e.g.  , 
by  a  combined  iteration  and  Newton-Raphson  method,  the  temperature 
T  and  the  mass  fractions  Y.  ,  i=  1,2,  3,5,  and  7  corresponding  thereto. 
Then  Eqs.  (II- 6.1 1)  are  solved  by  matrix  inversion  to  obtain  cr  ,  j  =  1 , 2,  3; 
the  integraL  F  in  Eq.  (II- 6.1 2)  can  then  be  evaluated  and  the  integration 
for  the  time  t  carried  out.  Note  that  when  Y  equals  the  equilibrium 
value  F~*oo  so  that  in  reality  this  final  value  does  not  have  to  be  known 
a  priori. 

In  Figure  II-5  there  is  presented  a  comparison  of  the  approximate 
and  exact  reaction  history  for  one  case;  shown  are  the  distributions 
of  Y3  ,  Yg  and  T  and  the  initial  conditions.  It  will  be  seen  that  in  general 
the  agreement  is  satisfactory.  As  a  result  of  these  considerations  it 
can  be  concluded  that  a  reasonable  understanding  of  the  interaction  of 
the  complex  mechanisms  involved  in  the  hydrogen-oxygen  reaction  is 
available;  of  course,  the  characterization  of  the  reaction  history  into 
two  time  periods  and  the  approximations  employed  in  the  simplified 
treatment  of  each  have  been  deduced  from  more  accurate  calculations. 

It  is  noted  that  the  extension  of  the  analysis  of  this  section  to  nozzle  flows 
according  to  either  the  streamline  or  one -dimensional  point  of  view  is 
straightforward. 
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FIG.  n-5.  Comparison  of  Approximate  and  Exact  Composition  and  Temperature 
Histories. 


7.  COMPARISON  WITH  EXPERIMENT  FOR  NOZZLE  EXPANSION 


As  mentioned  previously  there  are  available  limited  experi¬ 
mental  results  which  can  be  used  to  verify  the  overall  reaction 
mechanism  and  reaction  rates  employed  in  this  and  similar  studies 
of  the  hydrogen-oxygen  system.  In  addition  to  the  cited  experimental 
research  of  Schott  et  al.  ^  ^  there  are  available  the  tests  of 

nozzle  flows  of  the  combustion  products  of  hydrogen  and  air  carried 
out  by  Lezberg  and  Lancashire^  ^ ^ .  Of  particular  interest  herein 
are  the  measurements  of  the  temperature  along  the  axis  of  the  nozzLe. 
In  order  to  carry  out  a  comparison  between  experiment  and  theory  the 
pressure  distribution  along  the  axis  of  the  nozzle  had  to  be  estimated; 
the  approximate  transonic  theory  for  the  flow  in  the  throat  of  a  nozzle 
and  the  measured  pressure  distribution  along  the  nozzle  wall  were 
employed  to  obtain  the  estimated  pressure  distribution  shown  in 
Figure  II-6.  The  composition  in  the  settling  chamber  of  the  nozzle  was 
taken  to  be  close  to  equilibrium  at  the  specified  stagnation  pressure 
and  temperature  and  with  the  specified  equivalence  ratio.  The  equa¬ 
tions  of  Section  II-2,  with  the  reaction  mechanism  of  Eqs.  (II -3,1) 
with  RR  No.  2  selected  were  integrated  from  upstream  of  the  throat. 
The  composition,  velocity  and  temperature  distributions  along  the 
axis  were  computed,  Mr.  H,  Pergament  performed  these  calculations. 
Shown  on  Figure  II- 6  are  the  temperature  distributions  corresponding 
to  equilibrium  and  frozen  chemical  behavior  taken  from  reference  11-16 
and  corresponding  to  finite  rate  chemistry  computed  here.  Also  shown 
c.  r  the  two  experimentally  determined  temperatures  from  refer¬ 
ence  11-16,  In  view  of  the  approximations  attendant  upon  the  estimation 
of  the  pressure  distribution  and  in  view  of  the  possible  inaccuracies  in 
measurement,  the  agreement  with  respect  to  computed  and  measured 
temperatures  is  considered  satisfactory.  Clearly  additional  experi¬ 
mental  data  would  be  valuable  in  validating  the  reaction  mechanism 
and  rates. 
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PRESSURE  -  ATM 
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DISTANCE  FROM  NOZZLE  INLET -INCHES 

FIG.  II- 6.  Pressure  and  Temperature  Distribution  Along  the  Axis  of 


An  interesting  result  concerning  the  c or.c ent rat :ons  of  the  inter¬ 
mediates  in  a  non- equilibrium  flo  j  of  combustion  products  is  shown  in 
Figure  II-7a-c.  Shown  there  are  the  distributions  of  the  mass  fractions 
of  the  species  along  the  axis  of  the  nozzle  from  the  aforementioned 
calculation.  It  will  be  noted  that  recombination  of  H  and  O  taxes  place 
in  the  initial  part  of  the  expansion,  however,  downstream  of  the  point 
where  the  temperature  deviates  from  the  equilibrium  distribution  and 
where  the  mass  fraction  of  the  major  constituents  is  frozen.,  the  mass 
fractions  of  intermediates  continues  to  change;  of  particular  interest 
is  the  increase  in  H  and  the  descrease  in  OH.  This  is  somewhat 
unexpected;  in  the  extensive  studies  of  the  expansion  of  dissociated  air 
through  nozzles  and  around  bodies  in  hypersonic  flight  h  was  always 
found  that  the  concentrations  of  all  species  were  between  those  for 
complete  equilibrium  and  complete  frozen  flow.  The  behavior  here  tor 
the  hydrogen-air  system  appears  to  be  due  to  ;he  relatively  last  two- 
body  reactions  (1-4)  which  can  maintain  a  partial  equilibrium  of  the 
intermediates  while  the  major  constituents  are  essentially  frozen. 

This  indicates  that  perhaps  the  partial  equilibrium  analysis  discussed 
above  may  provide  reasonable  estimates  of  the  behavior  of  combustion 


products  during  nozzLe  expansion.  As  noted  above-  this  behavior  of 

(JT  6) 

H  and  OH  was  found  independently  by  Westenberg  and  Fac  in'  J 


8.  PRELIMINARY  CONCLUSIONS  OF  THE  ANALYSIS  FOR  CONSTANT 
PRESSURE  CHEMICAL  REACTION 

The  results  of  the  analysis  described  above  lor  constant  pres¬ 
sure  show  that  the  induction  time  decreases  rapid1  v  when  static 

_s 

temperature  and  pressure  increase.  Such  times  of  the  order  of  10 
sec  can  be  easily  obtained  for  high  Mach  number  flight  conditions  even 
at  high  altitudes.  As  a  consequence..  :r.  a  practical  combustor,  the 
ignition  zone  will  be  Limited  to  a  few  inches  after  mixing  and  the  energy 
release  zone  will.  be  of  the  same  order  fer  pressures  between  0  5  and 
1  atm.  On  this  basis  the  assumption  of  ir.viscid  flow  :or  the  combustion 
region,  where  large  variations  of  pressure  may  occur  represents  a 
good  approximation. 
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X  (DISTANCE  FROM  NOZZLE  INLET -INCHES) 

FIG.  II-7b.  Distribution  of  Mass  Fraction  Along  the  Axis  in  Nozzle  Flow. 
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For  such  conditions  the  analysis  of  the  diffusion  and  combustion 
process  along  a  given  streamline  can  be  divided  into  two  parts:  a 
mixing  region  where  transport  properties  are  important,  and  a  com¬ 
bustion  region  where  large  variation  of  state  occurs,  but  where  viscous 
effects  can  be  neglected.  In  general,  both  processes  occur  with 
pressure  variations. 

9o  APPLICATION  OF  INVISCID  ANALYSIS  TO  COMBUSTOR  DESIGN 

It  is  of  interest  to  consider  the  application  of  the  streamline 
anaLysis  discussed  above  to  actual  combustor  design.  The  objective 
herein  is  to  indicate  how  shock  waves,  which  could  be  associated  with 
heat  release,  can  be  avoided  by  proper  design  of  the  combustor,  i„e., 
by  expanding  the  flow  area  in  a  manner  compatible  with  the  progress 
of  combustion.  For  this  purpose  a  constant  pressure  fLow  will  be  dis¬ 
cussed  although  similar  considerations  apply  to  the  more  general  case 
of  pressure  varying  along  the  streamline. 

Consider  a  premixed  hydrogen-air  mixture  corresponding  to  an 
equivalence  ratio  of  unity,  with  a  constant  pressure  of  4,7  atmospheres, 
a  uniform  velocity  of  11,000  ft/sec..  .  and  an  initial  temperature  of  1240°K. 
Then  the  numerical  integration  for  the  reaction  history  as  discussed 
above  can  be  employed  to  compute  inter  alia  the  spatial  distribution  of 
mass  density  which  can  be  used  to  compute  the  distribution  of  the  area 
ratio  (A  p  )  required  to  maintain  constant  pressure  flow.  The  area 
ratio  A/Aq  where  Aq  is  an  arbitrary  flow  area  at  x=0  is  shown  in 
Figure  II-8;  it  will  be  seen  that  there  is  a  relatively  small  region  in 
which  the  density  decreases  rapidly  and  in  which  the  flow  area  must 
increase  if  constant  pressure  is  to  be  maintained.  The  total  increase 
in  cross  sec.ional  area  is  2,15  Note  that  if  a  different  analysis,  one 
based  on  constant  a  ea  were  employed-  a  significant  pressure  rise 
and  velocity  decrease  would  occur  in  the  region  of  rapid  heat  release; 
this  is  hown  e  g  .  in  reference  11-7. 
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FIG.  II- 8.  Distribution  of  Flow  Area  for  Constant  Pressure  Combustion. 


Now  consider  Figure  II - 9  which  indicates  how  this  result  can  be 
empLoyed  qualitatively  in  the  design  of  a  practical  two-dimensional 
combustor.  Assume  that  the  total  inlet  flow  area  is  divided  into 
streamtubes  and  that  the  mixing  between  hydrogen  and  air  is  controlled; 
in  particular,  assume  that  the  hydrogen  can  be  distributed  by  mixing 
non-uniformLy  across  the  inlet  to  the  combustor  as  would  be  the  case  if 
a  "swept"  injector  were  provided.  Now  select  the  length  AB  to  corre¬ 
spond  to  the  length  in  which  no  significant  change  in  area  ratio  is 
required;  at  B  turn  the  fLow  by  an  angle  a  so  that  at  C  the  original 
streamtube  has  increased  in  area  by  a  factor  of  2.15,  where  the  length 
AC  corresponds  to  the  Length  obtained  from  Figure  11-8,  beyond  which 
the  area  is  again  essentially  constant.  Note  that  a  2  1.15  (d/&  )  where 
d  is  the  original  streamtube  height  and  A  is  the  Length  over  wnich  the 
change  in  cross-sectional  area  is  effectively  required.  Now  for  the 
second  streamtube  assume  that  the  origin  for  the  chemicaL  history  is 
shifted  by  the  amount  BC  so  that  its  required  increase  in  cross- 
sectional  area  occurs  downstream  of  point  C  and  is  complete  in  the 
length  Ax  =  B  'C  '  =  BC,  If  this  process  is  carried  out  for  each  stream- 
tube,  there  wilL  be  obtained  the  duct  configuration  which  is  shown  in 
Figure  II- 9  and  which  invoLves  a  final  cross-sectional  area  2,15  times 
greater  than  the  original  one.  It  will  be  recognized  that  this  treatment 
corresponds  to  a  rough  streamtube  analysis  for  a  combustor.  This 
analysis  can  be  improved  and  refined  by  a  two-dimensional  treatment 
of  the  fLow  fieLd  with  finite  rate  chemistry.  An  analysis  thereof  is 
discussed  in  the  foLLowing  subsection. 

10.  METHOD  OF  CHARACTERISTICS  FOR  TWO-DIMENSIONAL 
FLOWS 

Consider  now  the  anaLysis  by  the  method  of  characteristics  of 
two-dimensional  supersonic  flows  with  finite  rate  chemistry.  It  will 


The  method  of  characteristics  for  non-equilibrium  gas  flows  corre¬ 
sponding  to  either  dissociation  or  vibrational  relaxation  has  been 
treated  in  references  11-20  to  11-23.  The  analysis  here  is  more 
general. 
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FIG.  II— 9 .  Schematic  Representation  of  a  C ortstant  Pressu re  C ombustion. 


be  shown  that  the  streamline  anaLysis  forms  ar.  integral  part  of  such  a 
method.  Although  for  practical  reasons,  two-dimensional  flow  is 
considered  explicitly  herein,  the  extension  to  the  axisymmetric  .case  is 
straightforward.  For  completeness  some  of  the  equations  of  Section  II-2 
are  repeated  here. 

The  equations  of  motion  in  coordinates  along  a  streamLine  s  and 
in  the  normal  direction  n  are 


(II- 10.1) 


a  s  on 


(II-10.2) 


The  equations  of  species  conservation  are 


9Y. 

pV  - -  =  w. 

P  3s  i 


(II- 10. 3) 


and  the  overaLl  continuity  equation  for  two-dimensional  flows  is 


i-  (pV)  +  pV  =  0  (II- 1 0.4) 

3  s  an 

These  equations  are  supplemented  by  the  requirements  of  energy  con¬ 
servation  [Eq.  (II-2.2)];  by  the  equation  of  state  [Eq.  (II- 2.5)];  by  the 
description  of  the  creation  terms  [Eqs.  (II-2.6)  and  (II-2.7)];  ar.d  finaLly 
by  the  enthalpy-temperature  relation  for  the  species. 

Now  foLLowing  the  usual  procedure  in  the  method  of  charac¬ 
teristics  (cf.  e.g.,  references  11-24  and  -25>i,  Eqs.  (II-2.2)  and  (II-2.5) 
are  differentiated  with  respect  to  s;  the  resulting  equations  and 
Eqs.  (Ef-10.1)  and  (11-10,4)  can  be  used  to  obtain 
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which  involves  onLy  the  partial  derivatives  9p/9s  and  90/8n  and  finite 
quantities.  As  usual  there  are  considered  the  total  differentials 


_1_  ^P  _  _1_  9p  ^  9p  ,  dn. 
p  ds  p  9s  p  9n  ds 


(II- 10.6) 


dQ  _  3_0  +  9^  .  dn, 
ds  9s  8n  ds 


(II- 10. 7) 


The  Eqs.  (11-10. 2)  and  (II- 1 0.5)  - (II- 1 0  ,7)  are  formally  considered  to 
yield  the  four  partial  derivatives  in  terms  of  finite  quantities  and  of 
d  In  p/ds,  d9/ds  and  dn/ds.  Select  dn/ds  so  that  the  determinant  of  these 
equations  vanishes;  there  are  obtained  two  values  of  dn/ds;  namely, 


dn 

ds 


R 


Wc 


— )  -1] 


2 


(H-10.8) 


The  condition  for  dn/ds  to  be  real  is  evident.  Moreover,  the  degener¬ 
ation  of  this  relation  to  the  special  case  usually  treated  in  gas  dynamics, 
i.e.  ,  constant  composition  and  constant  specific  heats  is  also  evident. 
FinalLy,  it  is  interesting  to  note  that  the  speed  of  sound,  which  arises 
quite  naturally  from  the  analysis,  is  associated  with  frozen  composition 
and  frozen  vibrational  energy;  namely. 


afs  =  (RT/W)[l-(R/cpW)]"  (II-10.9) 

Additional  restrictive  assumptions  must  be  made  a  priori  in  order  that 
there  enter  the  analysis  other  limiting  speeds  of  sound,  i.e.  ,  those 
associated  with  equilibrium  behavior  in  some  sense. 
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Denote  the  vaLues  of  dn/ds  corresponding  to  the  pLus  and  minus 
signs  by  (dn/ds)j  and  (dn/ds)^  ,  respectively;  introduce  and  pjj  where 

tan  Uj  =  {dn/ds) j 

(11-10.10) 

tan  Pjj  =  (dn/ds)^ 

The  vanishing  of  the  determinant  leading  to  the  special  values  of 
dn/ds  implies  along  the  line  (dn/ds)^  that 


where  d£  and  dn  are  differential  elements  along  each  characteristic 
line.  In  the  non  reactive  case  F  =  0  so  that  Eq.  (II  -  1 0.1 1)  and  (II- 1 0.1 2) 
reduce  to  those  usually  prevailing  in  gas  dynamics.  The  right-hand 
sides  of  Eqs.  (II -10.1 1 )  and  (11-10.12)  are  equivalent  to  the  terms  which 
arise  due  to  rotationality.  In  addition  to  the  two  characteristic  lines 
related  to  Uj  and  the  streamlines  are  also  characteristics  along 
which  the  reaction  equations  must  be  applied. 
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The  application  of  these  equations  to  a  unit  problem  within  the 
flow  field  as  shown  in  Figure  II -10  follows  the  chemical  numerical 
procedure.  Assume  that  at  two  adjacent  points  1  and  2  all  the  flow  and 
state  variables  are  known.  Then  from  Eqs.  (11-10*11)  and  (11-10.12)  in 
the  usual  way  two  finite  difference  equations  for  p3  and  03  can  be 
established  provided  the  finite  quantities  in  these  equations  are  evalu¬ 
ated  from  conditions  at  points  2  and  1  respectively.  Thus  the  location 
of,  and  the  values  of  p  and  9  at  point  3  can  be  found  to  a  first  approxi¬ 
mation.  Now  construct  point  4  on  the  streamline  inclined  at  03  passing 
through  point  3  and  at  94  passing  through  4  where  Linear  interpolation 
between  0a  and  9,  is  assumed;  further,  by  linear  interpolation  esti¬ 
mate  V,  p  and  Y;  ,  i=  1, 2  .  .  N-l  at  point  4.  Thus,  if  the  pressure  is 
assumed  to  vary  linearly  from  points  4  to  point  3,  these  estimates  pLus 
the  energy  equation  permits  the  conservation  Eqs.  (II -10.1)  and  (II-10.3), 
and  the  supplemental  equations  relative  to  the  state,  kinetics  and  species 
enthalpy-temperature  to  be  integrated  along  s  to  yield  the  composition 
at  point  3  in  the  same  fashion  as  in  the  usual  streamline  calculations 
described  above.  Iteration  can  also  be  used  if  required. 

It  is  pointed  out  that  the  procedure  outLined  here  assumes  that 
requirements  for  accuracy  in  integration  for  the  species  conservation 
equations  along  the  streamline  are  more  severe  than  the  requirements 
for  accuracy  in  the  flow  variables  p,  9,  V,  etc.  If  this  is  not  the  case, 
then  finite  difference  approximations  for  integration  along  the  stream¬ 
line  from  point  4  can  be  readily  employed. 

In  addition  to  the  unit  problem  for  the  flow  field,  special  unit 
problems  for  points  involving  shock  wave  and  for  boundary  points  must 
be  established.  For  the  flow  across  shock  fronts  use  of  the  von  Newman- 
Doring-Zeldovitch  jump  conditions  would  be  consistent  with  the  inviscid 
flow  assumption.  No  special  difficulty  with  respect  to  these  special  unit 
problems  would  appear  likely  in  principle  but  as  in  the  method  of  char¬ 
acteristics  for  classical  gas  dynamics,  difficulties  associated  with 
secondary  shocks,  numerical  accuracy,  programming  logic,  etc.  are 
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certain  to  arise  and  to  require  soLution.  An  analysis  employing  the 
method  of  characteristics  with  finite  rate  chemistry  along  the  lines 
described  here  will  be  valuable  in  combustor  and  exhaust  nozzle  design 
and  analysis. 

It  is  interesting  to  note  that  the  effects  of  viscosity,  diffusion 
and  heat  conduction  normal  to  the  streamlines  can  be  incorporated  in 
the  above  analysis  in  an  approximate  manner,  if  these  effects  are  slowly 
varying  with  respect  to  combustion  processes.  For  example,  in  the 
species  conservation  equations  the  effect  of  diffusion  normal  to  the 
streamlines  is  to  add  an  additional  term  9/9n(pD9Y  ,/9n) ;  an  average 
value  thereof  for  the  increment  along  the  streamline  from  point  4  to 
point  3  can  be  evaluated  in  finite  difference  form  from  a  knowledge  of 
the  state  at  points  1,  2,  and  4.  Similarly,  the  energy  equation,  which 
in  this  case  must  be  treated  in  differential  form,  will  have  an  additional 
term  arising  from  energy  transport  due  to  conduction  and  to  diffusion 
and  a  dissipation  term  due  to  viscosity;  again  a  finite  difference  repre¬ 
sentation  can  be  employed.  Finally,  the  viscosity  contribution  to  the 
streamwise  momentum  equation  can  be  taken  into  account.  For 
Laminar  flow  the  transport  coefficients  are  clearly  defined  whiLe  as  is 
well-known  for  turbulent  flow  effective  exchange  coefficients  must  be 
introduced. 


SECTION  III 

HETEROGENEOUS  TURBULENT  MIXING 

It  appears  that  under  flight  conditions  for  hypersonic  vehicles 
of  practical  sizes,  the  mixing  process  between  the  fuel  jet  and  the  air 
stream  will  be  predominantly  turbulent.  Accordingly,  the  present 
discussion  of  the  mixing  problems  connected  with  supersonic  diffusion 
fLames  will  be  exclusively  devoted  to  turbulent  flows.  It  will  be  realized 
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that  at  the  origin  of  mixing  there  will  always  exist  on  the  fuel  injector, 
boundary  layers,  one  in  air  and  one  in  the  hydrogen  fuel;  in  considering 
only  turbulent  processes,  it  is  assumed  that  either  one  of  these 
boundary  layers  is  turbulent  at  the  origin  of  mixing  or  that  transition 
from  laminar  to  turbulent  fLow  occurs  in  a  negligibly  short  distance 
downstream  from  the  origin  of  mixing. 

1.  GENERAL  REMARKS 

The  mixing  processes  in  supersonic  diffusion  flames  play  an 
important  role  in  determining  the  characteristics  of  a  combustor,  since 
the  fuel  and  oxidizer  must  be  brought  together  by  mixing  prior  to 
chemical  reaction.  As  mentioned  in  the  Introduction  the  relative  rates 
of  the  mixing  and  of  the  chemical  processes  taking  place  simultaneously 
in  the  fuel  jet  determine  the  characteristics  of  the  flow;  if  the  chemical 
reactions  are  relatively  fast,  i.e.,  for  the  hydrogen-air  system  of 
interest  herein,  if  the  initial  static  temperatures  and  the  static  pressures 
are  hig  i,  then  diffusion  will  be  rate  controlling.  In  a  limiting  sense 
under  these  conditions  the  behavior  corresponds  to  local  chemical 
equilibrium  and  heat  release  occurs  close  to  the  exit  of  the  fuel  injector. 
If  on  the  contrary,  the  chemical  reaction  is  relatively  slow,  the  mixing 
of  the  fueL  and  air  close  to  the  exit  will  correspond  to  a  heterogeneous, 
non-reactive  jet;  far  downstream  the  fuel  and  oxidizer  will  be  com¬ 
pletely  mixed  and  if  conditions  are  suitable  wilL  react  in  an  essentially 
inviscid  manner.  This  second  case  in  a  limiting  sense  corresponds  to 
"frozen"  mixing. 

The  two  limiting  cases  of  equilibrium  and  frozen  chemical 
behavior  are  the  easiest  to  treat  analytically  so  that  most  anaLyses  of 
viscous  reacting  flows  have  considered  only  such  cases.  However,  for 
combustor  design  it  is  necessary  to  determine  quantitatively  the  con¬ 
ditions  under  which  the  two  limiting  cases  prevail  and  the  characteristics 
of  flows  intermediate  to  the  limiting  cases,  i.e.  ,  of  fLows  with  finite 
rate  chemistry.  In  the  analysis  of  mixing  with  chemical  reactions  the 
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interaction  between  mixing,  combustion  and  boundary  conditions  intro¬ 
duces  severe  complications.  The  combustion  process  locaLLy  effects 
temperature  density  and  pressure  of  the  fLow.  If  the  flow  is  supersonic 
pressure  variations  propagate  in  the  flow  through  a  wave  mechanism, 
and  local  changes  of  static  pressure  can  influence  reaction  rates  and 
mixing  rates  elsewhere  in  the  flow.  These  coupling  effects  can  be 
important  in  practical  applications.  In  the  case  of  supersonic  diffusion 
flames  the  mixing  process  is  the  controlling  mechanism  and  a  strong 
increase  in  pressure  must  be  avoided  by  careful  selection  of  the 
boundary  conditions.  Therefore  it  appears  logical  to  analyze  in  the 
first  step  mixing  problems  where  the  pressure  can  be  assumed  to  be 
practically  constant  through  the  entire  flow  field. 

Even  in  the  case  of  a  constant  pressure  process  complications 
are  introduced  by  several  effects;  e.  g.  ,  by  the  details  of  the  flow  at 
the  origin  of  mixing;  in  the  boundary  layer  on  the  air  side  of  the  fuel 
injector  there  will  exist  regions  of  low  velocity  and  high  static  temper¬ 
ature  because  of  aerodynamic  heating.  Now  it  is  known  from  analyses 
such  as  those  presented  in  the  previous  section  that  the  reaction  history 
for  hydrogen-air  mixtures  is  sensitive  to  initial  temperature  as  well  as 
to  static  pressure.  Thus  the  chemical  times  and  lengths  associated 
with  the  static  temperatures  and  velocities  near  the  exit  of  the  fuel  jet 
may  actually  be  several  orders  of  magnitude  smaller  than  those  based 
on  the  static  temperatures  and  velocities  in  the  streams  external  to 
the  viscous  region.  If  this  is  the  case,  reaction  may  initiate  close  to 
the  fuel  jet  and  spread  rapidly  so  as  to  cause  extensive  amounts  of  heat 
release  in  flow  distances  shorter  than  those  expected  on  the  basis  of  the 
usual  overall  considerations  of  diffusion  and  chemical  times  for  idealized 
conditions . 

In  addition  to  the  effect  of  initial  boundary  layers  the  conditions 
under  which  limiting  chemical  behavior  can  be  expected,  are  com¬ 
plicated  by  the  wide  range  of  flow  velocities  which  can  occur  in  super¬ 
sonic  diffusion  flames  of  practical  interest.  For  example,  consider  the 
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case  wherein  the  fuel,  jet  has  a  small  velocity  compared  to  that  in  the 
external  stream.  In  this  case  chemical  reaction  may  first  appear  near 
the  axis  of  symmetry,  where  the  Lowest  velocities  occur,  and  may 
spread  rapidLy  outward  in  constrast  to  the  previous  discussion. 

The  problem  of  treating  analytically  the  simultaneous  effects 
of  turbulent  mixing  and  finite  rate  chemistry  is  formidable  if  some 
reaLism  is  desired.  The  preliminary  research  reported  here  considers 
first  the  turbuLent  mixing  of  heterogeneous  gas  streams.  Both  theo¬ 
retical  and  experimental  research  wiLL  be  described.  In  addition  some 
remarks  concerning  the  chemical  kinetics  of  turbuLent  reacting  fLows 
wiLL  be  made  to  indicate  an  essentiaL  difficulty  in  proceeding  to  the 
anaLysis  of  practically  interesting  cases  wherein  the  effects  of  finite 
rate  chemistry  are  important;  dearly  these  cases  are  the  intermediate 
ones  cited  above. 

The  problem  of  describing  quantitatively  turbuLent  shear  fLows 
has  occupied  the  attention  of  fLuid  mechanicians  for  many  years.  It  is 
generally  agreed  that  even  for  the  simplest  case  of  incompressible  flow 
the  available  knowledge  is  incomplete.  With  respect  to  turbuLent  shear 
fLows  involving  significant  gradients  of  density  the  present  state  of 
knowledge  is  poor.  In  view  of  this  situation  research  directed  toward 
establishing  some  of  the  engineering  fundamentals  of  turbuLent  mixing 
of  heterogeneous  gases  has  been  undertaken  at  PIBAL  and  at  GASL. 

The  flow  configuration  has  been  ideaLized  in  both  theory  and  experiment 
to  correspond  to  an  axisymmetric  fueL  injection  in  an  essentiaLLy  infinite 
oxidizing  stream.  In  a  practical  combustor  the  presence  of  walls,  of 
other  asymmetries  and  of  non- uniformities  in  the  oxidizing  stream,  e.g, 
due  to  gradients  in  stagnation  pressure,  wiLl  complicate  the  flow  so 
that  onLy  a  general  but  essential  understanding  of  the  phenomena  can  be 
expected  from  the  study  of  idealized  configurations  such  as  are  considered 
herein. 
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2,  THEORETICAL  ANALYSIS  -  GENERAL  EQUATIONS 

Consider  first  an  engineering  analysis  of  a  turbulent  jet  of  a  fuel 
such  as  hydrogen  in  an  air  stream;  assume  that  the  mean  turbulent 
quantities  are  described  by  the  equations  for  Laminar  flow  if  the  trans¬ 
port  coefficients  repLace  their  molecular  counterparts.  Then  (cf.  references 
III  -  1  red  -2)  for  constant  pressure  flows  the  following  conservation 
equations  are  obtained; 

Conservation  of  Momentum 


9u  9u  -l  9  ,  9u. 

pu  _  +  pv  _  -  r  — -  (£pr  —) 

dx  9r  9r  dr 


(m  -2.D 


Conservation  of  Energy 


9H  ,  8H  -l  9  r,  -l  9H> 
U  ax'  +  pv  9T  =  r  ¥7  [(DGpt  r  S7] 


.  i  n  -l  9u 

+  (Pt  -1)  P£Pt  ru  -gj 


(in-2.2) 


+  (St-1)  „,pt-  l  hiTfi] 
i=l 


Global  Conservation  of  Mass 


(pu)  +  r'1  ■—  (pvr)  =  0 


(in- 2.3) 


Species  Conservation 


r.  9Y .  „  9Y. 

-  +  pv  — -i-  =  r*  -f-  (peS  .r  — i-)  +  w  (111-2.4) 

;  or  9r  :  9r  1 


i=  1, 2.  -  N-L 


OI  j 

pu  -93T  +  pv 


3Yj  _  9 


i  a  .  9Y= 

r"  -57  (peSt’  r  ~ZT} 


(in-2.5) 


j=l,2.  -L 
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Eqs.  (III-2.1)  through  (III-2.5)  are  derived  from  the  general 
equations  of  turbulent  motion  by  introducing  the  concept  of  eddy 
viscosities  that  permits  to  determine  the  contribution  of  the  fluc¬ 
tuating  quantities  of  the  flow  as  a  function  of  the  local  gradients  of  the 
average  properties.  On  this  basis  the  transport  properties  due  to  the 
turbulent  motion  can  ho  expressed  in  the  same  form  as  the  transport 
properties  due  to  the  molecular  motion.  However,  the  coefficients  of 
the  terms  that  represent  the  molecular  motion  are  functions  only  of 
local  physical  properties  of  the  gas,  while  the  coefficients  related  to 
the  turbulent  motion  usually  are  related  to  gradients  of  physical 
properties  of  the  gas.  When  the  contribution  due  to  molecular  motion 
is  small  and  can  be  neglected,  then  the  eddy  viscosities  are  the  same 
for  all  the  transport  properties;  this  corresponds  to  the  turbulent 
Prandtl  and  Schmidt  number  equal  to  unity.  In  the  case  that  the  con¬ 
tribution  of  the  molecular  transport  is  important  and  must  be  retained, 
then  it  is  represented  in  the  equations  in  an  average  form,  by  intro¬ 
ducing  correcting  factors  to  the  different  transport  properties.  This 
is  done  by  assuming  that  the  turbulent  Prandtl  and  Schmidt  numbers 
are  different  from  unity  but  constant  at  any  given  cross  section  of  the 
flow.  On  the  equations  presented,  a  single  turbulent  diffusion  coef- 
ficient  D  corresponding  to  a  single  turbulent  Schmidt  number  St  =  jj 
has  been  used. 

Fundamental  to  the  application  of  these  equations  to  the 
problem  of  a  turbulent  jet  is  a  requirement  to  describe,  at  least 
approximately  the  compressible  eddy  viscosity  e.  In  the  past  theo¬ 
retical  analyses  of  turbulent  compressible  mixing  (see  references  HI-3 
and  -  4)  have  carried  over  the  Prandtl  hypothesis  generally  accepted 
for  incompressible  flow;  namely,  that  the  eddy  viscosity  is  essentially 
uniform  across  the  mixing  region  and  varies  at  most  with  the  stream- 
wise  coordinate. 

A  different  analysis  for  compressible  mixing  problems  has  been 
presented  by  the  second  author  in  reference  IH-1.  In  the  analytical 
work  presented  here  the  approach  presented  in  reference  ni-1  will  be 
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r3  = 


cu 


and 


£  r3  r3  u 

p.3u.p 

J  J 


ep  -  (ep). 


(III- 3. 5) 


(III -3. 6) 


2  2 

Assume  (ep)  =  f(x),  then  near  the  axis  =  f(x)  and 


Eq.  (Ill- 3. 3)  becomes 


W 


if  U=u/u.  and 
J 


3u  =  ^  ^c  _9  ~  9u 

9x  ii'  p^  9iJ  9  5 


{£p)c 

p.u.a 
J  J 


dx 


then 


9U  _  a  _9_  ,r  9 U. 
9?  5  9?  U  95 


(III- 3.7) 


(in- 3. 8) 


(IH-3.9) 


The  expression  (IH-3.6)  is  accurate  near  and  far  away  from  the 
axis,  and  is  only  an  approximate  expression  moving  away  from  the  axis. 
The  error  introduced  is  due  to  the  fact  that  terms  of  the  order  (~~4~-)  5  3 
are  neglected  with  respect  to  one  in  Eq.  (III-3.5).  Indeed 


9t|i 


pu  =  (pu)  [1+  —  - — — 
c  4  ( pu) 


9s  (pu), 
9*a 


5s  + 


(III- 3. 10) 


Then,  if  the  second  term  is  neglected  in  Eq.  (Ill -3.5) ,  Eq.  (HI-3. 6)  gives 


(III-3.il) 


~  2  ~ 

epaurs  ~  (e p)  <!>  =  f(x)  iji  s 

This  expression  Leads  to  Eq.  (HI-3. 9).  The  soLution  of  Eq.  (III-3.9)  has 
been  indicated  in  references  IH-1  and  -5  ,  and  is  independent  of  the 
energy  equation. 


4.  EVALUATION  OF  THE  EDDY  VISCOSITY 

In  reference  IH-1  the  suggestion  has  been  made  that  the  expres¬ 
sion 


p  e  =  k  (u  -u  .)p(r,). 
c  max  min  c  4  i 


be  used  for  definition  of  the  eddy  viscosity  e,  where  (rO;  is  the  width  of 

2 

the  mixing  zone  in  the  incompressible  plane  wherein  the  velocity  changes 

from  u  to  (u  +  u  )  A  .  The  value  suggested  for  the  constant  k  was 
max  max  c  2  66 

0.025  as  given  by  PrandtL,  then 


p  e  =  0.025  u  1 1  -u 
pc  e  c 


(n).  p 
'  fine  e 


(m-4.1) 


This  expression  indicates  that  the  turbulence  would  decay  when  the  two 
jets  have  the  same  velocities  even  if  the  two  streams  have  different 
temperature  or  composition.  This  is  not  consistent  with  the  fact  that 
in  these  cases  a  dissipative  mechanism  stilL  exists  due  to  heat  conduc¬ 
tion  or  concentration  changes  that  would  sustain  turbulence.  It  would 

aLso  indicate  that  a  sharp  decrease  of  mixing  rates  would  occur  for 
u; 

—  1  with  rapid  changes  of  mixing  rates  in  both  sides  of  the  con- 
ue  uj  ... 

ditions  — ■*-  ^  1  .  A  modification  of  the  expression  for  the  eddy  viscosity 
ue 

has  been  suggested  by  the  first  author  in  reference  Ill-b.  The  expres¬ 
sion  used  here  has  the  form 


pe  =  (pe)c  =  kr^  lpeue-pcuc  I 
z 


(IH-4.2) 
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where  k  is  the  constant  given  by  PrandtL  and  is  assumed  to  be  equal  to 
0.025.  This  expression  still  follows  the  basic  concept  of  mixing  length 
of  Prandtl  and  the  basic  concept  of  reference  III  - 1  with  respect  to  the 
dependence  of  eddy  viscosity  on  spacial  coordinates.  In  this  expression 
the  quantity  r^  is  the  width  of  the  mixing  zone  in  the  physical  plane, 
defined  as  the  value  where  the  quantity  pu  has  an  average  value  between 
the  value  at  the  axis  and  the  vaLue  in  the  outside  flow.  Therefore 

*  ,Pc“c't'e“e>  I 

2 

It  has  been  found  that  this  expression  is  in  better  agreement  with  the 
experimental  resuLts  presented  here. 

5.  THE  FREE  JET  IN  A  QUIESCENT  ATMOSPHERE 

Consider  next  the  problem  of  free  jet  in  a  quiescent  atmosphere. 
This  problem  is  of  interest  because  it  has  been  investigated  experi¬ 
mentally  for  a  large  variety  of  cases.  Assume  that  the  boundary  layer 
at  the  origin  of  mixing  (x=  §  =  0),  is  negligible.  Then  U(?  ,  iji )  is  to 
satisfy  the  conditions 


U(0 ,  i|i)  =  1  0  <  ij/  <  a 

=  u  /u.  iji  >  a  (III- 5.1) 

lim 

7  U(5  1 1)  =  u  /u. 

-  co  y  e'  j 

The  solution  of  the  problem  specified  by  Eqs.  (III-3.il)  and  (III -5.1)  is 
a  standard  one  in  the  theory  of  heat  conduction  (cf.  references  HJ-7 
and  -8).  The  solution  is  given  by  the  equation 
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[U-(u  /a.)]  x 

u  =  TTT-Vt!  =  (2?/a)1  exp[4/a)s/4(?/a)]  f  exp[ -  r  '  a/4(e/a) ] 

[ 1  -  (ue/uj  ]  o 

(IH-5.2) 

IQ[r  '  it /a.)  Ue/a.)  -1  ]r  '  dr  ' 

and  has  been  tabulated  by  Masters f111'7)  in  terms  of  the  function  P* 
where  U  =  P*(Z/v  ,  R/y  ) Uc  with  Z/y  =  (2?/a)  R/y  =  (*/a)  (Z/v) .  Of 

particular  interest  in  the  present  discussion  is  the  velocity  decay  along 
the  axis  in  terms  of  5  ;  there  is  obtained  from  the  solution  that 


Uc  =  1  -  cxp(-a/4e) 


(HI-5.3) 


For  the  case  u  =0, 
c 


(ep)  -  f(x)  =  kp  u  r , 

C  C  C  I 

2 


(HI -5. 4) 


(pu>r,  =  2  Pcuc 
2 


(III- 5. 5) 


Then  the  quantity  is  defined  by 


r  i  2p,u 

(r  >3  =  r  — J— 

2  L  o  pu 


(IH-5.6) 


Asymptotically  (£  —  00)  this  expression  becomes  equal  to 


r  1  — M-^-)  ? 
2  °e 


(HI-5.7) 
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The  quantity  asymptotically  becomes 


u  .a 


and 

P  1 

[ ( p e)  ]  =  k  p  u.a  (— )  =  constant 

c  as  e  J  Pe 


(III- 5. 8) 


(III-5.9) 


Mr.  G.  Kleinstein  has  found  from  the  analysis  of  numerical  and 
experimental  data  that  the  quantity  (ep)^  remains  constant  in  a  Large 
zone  of  the  mixing,  and  that  the  asymptotic  expression  can  be  used  with 
good  approximation  also  in  the  region  of  the  fLow  near  the  exit  of  the 
jet.  On  the  basis  of  this  observation,  the  evaluation  of  Eq.  (Ill  —  3 . 9) 
becomes  very  simple.  If  the  quantity  k  is  assumed  to  be  the  constant 
given  by  Prandtl,  by  introducing  Eq.  (III-5.9)  in  Eq.  (HI-3, 10)  it  results 


—  =  40  (— j-)  2  +  constant 

a  n  a 

Me 


(III- 5. 10) 


x. 

Denote  the  constant  in  this  equation  as  — ,  then  Eq.  (III-5.3)  becomes 

cL 

ue  =  1  -  exp  {-10[~  *  (in-5.11) 

e  e 

where  a^  is  an  effective  area  of  the  jet. 

P.  I 

a  =  a(-J-)  (HI-5.12) 

e  0 

e 

This  resuLt  is  in  accord  with  the  observation  of  Thring  and  Newby 

that  a  homogeneous  compressible  jet  has  the  same  behavior  as  an 

incompressible  one  provided  the  actual  jet  radius  a  is  replaced  by  an 

effective  radius  a  . 

e 
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The  accuracy  of  this  simple  result  can  be  assessed  as  folLows: 
Consider  experiments  with  u  =0  yielding  the  decay  of  center  line 
velocity  (x/a);  compute  from  Eq,  (III -5.3)  for  each  u^  the  value 

of  e/a  and  compare  the  computed  vaLues  with  the  measured  values 

i 

plotted  in  terms  of  0„025(P  /p  ,)2(x/a).  This  has  been  done  in 

^  J 

Figure  III—  1 ;  there  have  been  considered  data  for  homogeneous, 
isothermal  (i  e,,  incompressible)  jets  (references  III  -  9  r  -10,  -11, 
and  -13),  isothermal  jets  of  carbon  dioxide  (reference  III- 13),  heated 
homogeneous  jets  (reference  III— 11),  and  high  speed  homogeneous  jets 
(reference  III  -10).  A  density  range  of  A  .  p  fo  ;<2is  represented  by  the  data. 
With  the  singLe  exception  of  the  high  Mach  number  data  of  Warren 
(Mj  =  2.60)  the  correlation  of  the  results  is  as  good  as  can  be  expected 
It  is  noted  that  a  similar  discrepancy  with  respect  to  these  data  was 
found  in  reference  III  - 1  and  remains  unexplained,  but  that  for  M.=  1.51 
the  agreement  is  satisfactory  herein. 

The  results  of  Figure  III- 1  imply  that  the  constant  in  Eq,  (III-5.10) 

is  approximately  -0.175  independently  of  the  ratio  p  /p.  ;  thus  the  trans- 

~  ^  J 

formation  from  the  5  ,  Jj  pLane  to  the  x,  it  plane  is  for  this  case  of  ue  =  0 
simply 

2  =  40(fl)2  [1  +  0.175]  (HI-5.13) 

a  °e  a 


Note  that  the  final  transformation  from  x,  (  to  x,  r  involves  application 
of  Eq.  (Ill - 3  4)  and  requires  specification  of  the  density  distribution  at 
least  in  terms  of  S  and  -j  , 

In  Figures  HI-2,  -3,  -4,  and  -5,  some  comparisons  between 
measured  and  calculated  quantities  for  the  case  of  heterogeneous  jets, 
investigated  experimentally  in  reference  III -13,  are  shown. 

Figures  IH-2,  -3,  and  -4  give  the  decay  of  velocity  and  concen¬ 
tration  along  the  axis  for  nitrogen,  carbon  dioxide,  and  helium  jets. 

In  the  case  of  helium  jet  the  experimental  velocity  distribution  shown 
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X  =  X/r 

Q 


FIG.  Ill- 2.  Comparison  of  Theoretical  and  Experimental 
Distributions  of  Center  Line  Velocity  and 
Concentration  for  Nitrogen  Jet. 
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X  =  X/rQ 

FIG.  Ill- 3  .  Comparison  of  Theoretical  and  Experimental 
Distributions  of  Center  Line  Velocity  and 
Concentration  for  Carbon  Dioxide  Jet. 


ref  nr -13 


X  =  X/rQ 


FIG.  Ill  -  4 .  Comparison  of  Theoretical  and  Experimental 
Distributions  of  Center  L.ine  Velocity  and 
Concentration  for  Helium  Jet, 
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i.O 


FIG.  Ill  -  5 .  Comparison  of  Theoretical  and  Experimental  Radial 
Distributions  of  Concentration  for  Helium  Jet. 


does  not  agree  with  the  analytical  value'  however,  it  is  stated  in  refer¬ 
ence  III- 13  that  "velocity  data  could  not  be  obtained  for  the  jet  of  heliurn 
because  of  the  low  impact  pressures  that  it  produced." 

In  Figure  III-5,  a  radial  profile  distribution  for  the  helium 
experiment  is  shown.  The  agreement  with  the  analysis  is  good. 

6.  ANALYSIS  OF  MIXING  FOR  P  AND  St  DIFFERENT  FROM  UNITY 

In  order  to  determine  concentration  and  temperature  distribution 
in  the  flow,  Eqs.  (III-2.4),  (III-2.5),  and  (III-2.2)  must  be  solved.  These 
equations  involve  the  turbulent  Prandtl  and  Schmidt  numbers. 

If  these  numbers  are  unity,  then  for  jet  flows,  generalized  Crocco 
relations  exist  among  the  velocity  stagnation  enthalpy  and  eLement  mass 
fractions  so  that  the  velocity  distribution  determines  the  entire  flow.  If 
the  flow  is  nonreacting  Eqs.  (III-2.4)  and  (IH-2.5)  are  Identical,  there¬ 
fore,  the  same  relations  are  valid  also  for  species  mass  fractions. 

In  the  case  of  jets  discharging  in  quiescent  atmosphere  is  zero 
and  if  Pt  and  St  are  unity  the  ratio  ~  and  Kcj  (concentration  of  the 
constituent  of  the  jet)  are  exactly  the^same. 

For  the  case  when  and  are  different  from  unity,  solutions  of 
Eqs.  (Ill- 2.4)  and  (HI-2. 5)  can  still  be  carried  over  directly  to  obtain 
concentration  and  temperature  provided  P  and  St  are  constant  or  at  most 
functions  of  the  streamwise  coordinate. 

Consider  a  solution  to  Eq.  (HI-2.4)  applied  to  one  foreign  species, 
i.e.,  for  Y.sK  and  for  w,  =0,  i.e.,  for  no  chemical  reaction.  Application 
of  the  von  Mises  transformation  and  the  previous  considerations  concerning 
e  leads  with  no  further  approximation  to 


3K 


(ec) 


cfx  ~  (~ 


(III- 6.1) 


69 


Now  the  available  data  on  turbulent  Schmidt  numbers,  ,  for  a  variety 
of  mixtures  indicates  that  an  assumption  of  a  constant  value  is  therefore 
i  ason?b'e  (cf.  references  III- 1  2  and  -14).  Accordingly,  for  these 
jet  flows  assume  to  be  constant  and  introduce  a  new  independent 
variable  £ ^  such  thatv 


e 

-l 


(ep) 


p.u. 

J 


dx / 
a 


Then  Eq,  (III —6*1)  becomes 


(IH- 6.2) 


JK  =a  a^jBK) 
3?  i[i  3i|(  3  i!t 


(III- 6. 3) 


Eq.  (Ill -6.3)  has  the  same  form  as  Eq.  (III-3.il).  For  a  jet  discharging 
in  quiescent  atmosphere 


K(0,  y)  =  1  0  <  \li  <  a 

=  0  J  >  a  (01-6.4) 

_Um  K(?,*)=0 
•j  -  oo 

Then  the  solution  for  K  is  the  same  as  that  for  U  given  by  Eqs.  (Ill  —  5.2) 
and  (III- 5, 3)  provided  5  is  replaced  by  1  £  ;  clearly  if  Z  1 ,  K  ==  U, 

i.e.,  a  Crocco  relation  exists  between  velocity  and  concentration  of 
foreign  gas.  The  concentration  along  the  axis  is  from  the  solution  given 
by 


Note  that  with  no  difficulty  S 


could  be  taken  to  be  a  function  of  x. 
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Kc  =  1  -  exp[-St/4(?/a)] 


(III- 6.5) 


while  from  reference  HI-8  the  concentration  profiles  are  given  as 
K/Kc  =  P*(Z  '/y,R'/y),  where  Z  '/y  =  (2^  /a) R/y  =  (*/a)Z  '/y. 

For  isothermal,  isobaric  conditions  pW  x=  constant  or  more 
precisely 

d/p ;  =  fK+  (W ./w  )(1  -K)]"1  (in-6.6) 

Thus  from  Eqs.  (III-6.5)  and  (III-6.6)  evaluated  at  the  center  line, 

Eq.  (III-3.5)  can  be  used  to  carry  out  the  transformation  back  to  the 
x  -  it  pLane.  Eq.  (Ill- 6.2)  for  ug  =  0  and  the  approximate  reLation 

P.  1 

(pe) c  =  0.025  u.a  Pe  (-^-) 

gives 

0  i 

?=  S'  [0.025  (-£)  x- constant]  (III- 6.7) 

t  p. 


Therefore,  from  experimental  data  of  Kc=Kc(— )  it  is  possible  to  deter¬ 
mine  experimentally  the  value  of  .  The  results  of  these  considerations, 

applied  to  the  jets  of  helium  and  carbon  dioxide  of  Figures  III-3  and  -4, 

Pe 

indicate  that  in  this  case  S  ~  1.  Note  that  a  change  in  density  ratio  — 

t  6  7  P- 

of  over  an  order  of  magnitude  is  represented  here.  J 


A  simiLar  analysis  can  be  performed  when  the  temperature  of  the 
jet  is  different  from  the  temperature  of  the  surrounding  flow.  In  this 
case 


p/o.  -  T ./T 
'  J  J 


so  that  the  temperature  field  must  be  computed.  For  a  constant 
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turbulent  Prandtl  number,  not  unity  in  general,  the  previous  analysis 
for  K  can  be  carried  through  to  yield  the  temperature  field;  the 
dependent  variable  becomes 


=  -  (T  -  Te)(T  -Te)_1  (III -6.9) 

while  the  independent  variable  becomes 

?a  5  i'*  (p_/ pJ ("e/u.a)  dx'  =  «  ?  (HI-6.10) 

1  “0  e  J  J  L 

The  solution  with  ft  =  1  at  x=?2  =  0  is  again  given  by  Eqs.  (Ill -5. 2)  and 
(IH- 5. 3)  provided  is  replaced  by'  F  and  U  by  P. 

7.  EXPERIMENTS  CONCERNING  HETEROGENEOUS  MIXING  OF  TWO 
MOVING  STREAMS 

In  order  to  obtain  fundamental  data  on  the  mixing  processes 
involved  when  there  is  an  external  flow,  a  series  of  experiments  con¬ 
cerning  the  axisymmetric  mixing  of  hydrogen  and  air  have  been  carried 
out  by  Mr.  L.  Alpinieri  at  PIBAL  it  a  mixing  rig  utilizing  the  super¬ 
sonic  facility  at  PIBAL;  the  static  temperatures  were  roughly  ambient 
so  that  no  chemical  reaction  occurred.  Two  s  ,ts  of  experiments  were 
performed;  in  one  a  subsonic  hydrogen  jet  exhausted  into  an  axisym¬ 
metric  stream  of  air  with  a  Mach  number  of  3.0,  The  experimental 
setup  is  shown  in  Figure  III-6.  In  the  second  a  subsonic  jet  of  hydrogen 
exhausted  into  a  subsonic  air  stream;  the  experimental  rig  is  shown  in 
Figure  III-7  In  both  sets  of  tests  the  distributions  of  static  and  total 
pressures  and  of  hydrogen  concentration  throughout  the  mixing  region 
were  measured-  The  pressures  were  measured  by  standard  techniques 
while  the  hydrogen  concentration  was  measured  by  a  thermal  con¬ 
ductivity  cell  which  was  calibrated  against  Known  hydrogen-air  mixtures. 
In  the  sub  sonic  -  super  sonic  tests  the  stagnation  pressure  of  the  air 
stream  was  1.3  atmospheres-  In  both  sets  of  tests  the  stagnation 
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temperatures  of  all  streams  was  approximately  290°K.  The  experi¬ 
mental  results  and  comparisons  with  theory  will  be  given  later  in 
extenso;  for  the  present  Only  the  results  indicating  some  important 
points  will  be  presented,  i.  e.  ,  those  pertaining  to  the  importance  of 
initial  boundary  layers,  to  isovelocity  mixing,  and  pressure  gradients. 


For  this  purpose  consider  first  the  distributions  of  concen¬ 
tration  obtained  in  the  subsonic-supersonic  tests;  shown  in  Figure  III-8 
is  the  distribution  of  mass  fraction  of  hydrogen  denoted  K  along  the 

center  line  for  four  different  values  of  the  velocity  ratio,  u./u  =  U.  . 

J  e  J 

Shown  in  Figure  III-9  are  the  radial  distributions  of  K  for  the  same 

velocity  ratios  and  for  two  downstream  stations  x/a  =  13  and  32.4.  It 

will  be  noted  first  from  these  figures  that  the  mixing  of  hydrogen  and 

air  is  rapid,  i.  c.  ,  drops  to  0.1  with  a  distance  corresponding  to 

x/a<  12  for  u.  =  0.36ue.  These  results  are  qualitatively  in  agreement 

with  the  anaLysis  discussed  above.  For  a  given  u  ,  as  u.  increases 

b  e’  j 

from  zero  the  theory  predicts,  at  a  given  station  ^  ,  an  increase  in  K^. 

For  the  form  of  the  eddy  viscosity  proposed  here  the  increase  of  Kc  is 

gradual  in  the  region  of  u.  =  u  ,  however  no  changes  in  sign  of  the 

J  e 

variation  of  Kc  with  uj  is  predicted.  This  is  not  true  for  the  form  of  the 
eddy  viscosity  proposed  in  reference  III-l,  which  would  indicate  an 
increase  in  Kc  until  uj  =  ue  ,  and  then  a  decrease  of  Kc  for  values  of 

Uj>Ue. 


A  direct  comparison  between  these  experiments  and  anaLysis 
requires  the  development  of  analyses  that  takes  into  account  the 
presence  of  the  boundary  layer  on  the  air  side  at  the  origin  of  the  mixing. 
These  analyses  can  be  derived  approximately  by  extending  the  present 
analysis  if  the  thickness  of  the  boundary  layer  is  not  too  large.  The 
presence  of  the  air  boundary  Layer  for  the  experiments  involving  mixing 
of  air  and  hydrogen  is  very  important,  because  the  momentum  flux  in 
the  air  boundary  Layer  is  very  large  with  respect  to  the  momentum  flux 
of  the  hydrogen  stream,  due  to  the  low  densit.v  of  the  hydrogen.  In  these 
experiments  for  the  data  presented  here,  the  boundary  layer  is  very 
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FIG.  m-9.  Radial  Concentration  Profiles  for  Various  Velocity  Ratios; 
Subsonic  -  Supersonic  -  Mixing  x/a  =  32.4 
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large,  the  thickness  is  of  the  order  of  one-half  inch  ,  therefore  the 
momentum  flux  deficit  in  the  boundary  layer  is  Larger  than  the 
momentum  of  the  hydrogen.  It  is  difficult  to  reach  even  qualitative 
conclusions  from  these  experiments-  and  the  resuLts  shown  in 
Figures  III-8  and  -9  should  not  be  interpreted  in  terms  of  step 
distributions  of  velocity  at  the  origin  of  mixing.  Indeed,  it  is  of 
interest  to  note  the  following:  The  great  density  difference  between 
hydrogen  and  air  impLies  that  the  boundary  layer  in  the  air  stream  at 
the  origin  of  mixing  will  generally  play  an  important  role  in  the  mixing 
processes  as  well,  as  mentioned  previously,  in  the  chemical  kinetics 
under  practical  conditions  of  interest  in  supersonic  diffusion  flames. 

One  of  the  main  reasons  for  setting  up  the  rig  for  subsonic  - 
subsonic  mixing  was  to  reduce  significantly  the  initial  boundary  layers. 
Indeed,  the  results  of  subsonic-subsonic  mixing  between  hydrogen  and 
air  indicate  that  the  qualitative  agreement  between  theory  and  experi¬ 
ment  is  not  due  to  the  effects  of  the  initial  boundary  Layer.  In  this  case 
the  boundary  Layer  in  the  air  stream  at  the  origin  of  mixing  is  thin; 
consequently,  the  effect  of  the  initial  boundary  layers  is  negligible  at  a 
small  distance  downstream  of  the  origin  of  the  mixing  region, 

A  comparison  between  experimental  and  analytical  resuLts  for 
axial  distributions  of  concentration  of  hydrogen  for  three  ratios  of  jet 
to  free  stream  velocities  is  presented  in  Figure  III-10.  In  Figure  III -11 
a  comparison  of  radial  distributions  for  a  velocity  ratio 
shown  for  values  of  x/R  =  5  and  10.5  „ 

The  trend  given  by  the  analysis  is  in  agreement  with  experi¬ 
ments.  The  concentration  at  any  given  station  x  =  constant  increases 
graduaLLy  when  uj  changes  from  0.7  to  1-3  when  the  mixing  extends  to 
the  axis.  The  concentration  initially  changes  very  rapidly  with  distance 
and  then  when  smalL  vaLues  of  concentration  are  reached  changes  slowly 
with  distance.  Both  resuLts  indicate  that  the  mixing  process  is  very 
rapid  also  for  vaLues  of  ~1  and  the  shape  of  given  by  the  analysis 
appear  to  be  in  agreement  with  experiments. 


,./u  =  1  is 
J  e 
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FIG.  Ill- 11.  Comparison  Between  Experimental  and  Analytical  Results 
for  Radial  Distributions  of  Concentration  of  Hydrogen 
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8.  EFFECT  OF  TURBULENCE  ON  REACTION  RATES 


In  order  to  establish  properties  of  mixing,  reacting  flows;  it  is 
necessary  to  combine  mixing  analyses  of  the  type  described  herein  with 
descriptions  of  the  chemical  kinetics  in  turbulent  flows.  Formally, 
expressions  for  the  creation  terms  w.  in  the  species  conservation 
equations  are  required.  It  will  be  recalled  that  the  equations  employed 
here  to  investigate  mixing  are  in  terms  of  time-averaged  quantities  so 
that  the  desired  descriptions  are  of  the  time -averaged  mass  rate  of 
production  of  each  species  per  unit  volume  per  unit  time.  While  the 
effect  of  turbulence  in  increasing  ihe  mixing  of  fuel  and  oxidizer  is 
understood  from  a  phenomenological  point  of  view,  the  effect  of  turbu¬ 
lence  on  kinetics  is  not  clear. 


The  effect  of  turbulence  on  the  rate  of  flame  propagation  has  been 
extensively  studied  but  apparently  the  fundamentals  of  the  effect  are 
stilL  subject  to  debate.  In  reference  III-15,  Hawthorne,  Weddell,  and 
Hottel  considered  the  effect  of  turbulence  to  be  described  by  an 
"unmixed  factor"  which  indicates  the  degree  of  inhomogeniety  of  fuel 
and  oxidizer;  they  considered,  however,  that  the  actual  kinetic  rates  of 
chemical  behavior  were  infinite  so  that  fueL,  oxidizer  and  products  were 
in  equilibrium  locally  with  respect  to  space  and  time.  Scientific  investi¬ 
gation  of  the  chemical  kinetics  of  turbulent  flows  according  to  statistical 

,  ,  ,  .  .  .  ,  t  ,  .  ...  lrr-  .  (Ill- 16,  -17)] 

theory  has  been  initiated  by  turbulence  specialists  |cf.  borrsin  J 

but  the  idealizations  necessary  for  progress  are  so  drastic  as  to  provide 

little  hope  for  application  to  chemical  systems  of  practical  interest  in  the 

near  future. 


It  is  the  purpose  of  this  discussion  to  indicate  qualitatively  the 
influence  of  turbulence  on  chemical  kinetics.  For  this  purpose  it  is  useful 
to  recall  the  basis  for  the  rate  constants  employed  in  the  analysis  of 
Section  II  and  shown  explicitly  in  connection  with  Eqs.  (II  -  3  1).  These 
rate  constants  have  been  obtained  by  theoretical  and/or  experimental 
means  under  the  cited  assumption  that  thermodynamic  equilibrium  roughly 
prevails  and  thus  that  '.on-equilibrium  prevails  only  with  respect  to  chemical 


81 


composition.  Therefore,  it  is  assumed  that  equipartition  of  energy 
among  particLes  and  internal  degrees  of  freedom  exists  at  Least 
approximately.  Now  in  high  speed  turbulent  flows  these  assumptions 
may  be  open  to  question;  therefore,  substituting  classical  rate  expres¬ 
sions  into  the  creation  terms  interpreted  as  instantaneous  rates  and 
the  time  averaging  thereof  in  order  to  obtain  time -averaged  values  for 
vc  may  not  be  vaLid.  However,  it  appears  that  the  following  statements 
can  be  made:  If  at  a  fixed  point  in  the  fLow,  the  absoLute  velocity 
squared  fluctuations  which  may  arise  from  either  velocity  or  temper¬ 
ature  fluctuations  are  significant  with  respect  to  2kT/mj  ,  then  these 
fluctuations  influence  significantly  the  collision  processes.  This 
appears  to  be  the  case  for  mixing  processes  of  very  high  velocity 
streams  having  very  large  ratio  of  static  to  total  enthalpy.  In  this  case 
the  enthalpy  connected  with  the  fluctuating  process  is  of  the  same  order 
of  the  static  enthalpy.  Therefore  the  mechanism  of  collisions  and  the 
number  of  a  given  type  of  collision  can  be  quite  different  from  the  vaLues 
given  by  statistical  mechanics  on  the  basis  of  thermodynamic  equilib¬ 
rium.  The  effect  of  the  turbulence  on  the  mechanism  of  collision  can 
be  discussed  qualitatively  by  considering  a  reaction  step  following  an 
Arrhenius  law.  The  turbulence  effects  the  value  of  T  on  the  exponent 
of  the  expression  and  the  function  of  T  which  multiply  the  exponential 
function.  Consider,  for  example,  a  reaction  step  following  an 
Arrhenius  law  [reactions  1-4  to  the  right  or  to  the  Left  and  reactions 
5-8  to  the  left  in  Eqs.  (II -3.1)  ] .  The  reaction  rate  in  this  case 
increases  greatly  with  the  reLative  velocity  between  colliding  particLes 
as  evidenced  by  the  exponential  dependence  of  the  rate  constant  on 
temperature.  Therefore,  it  can  be  expected  that  a  large  fluctuation 
in  the  absolute  velocity  squared  will  increase  the  reaction  rate  for 
reaction  steps  following  an  Arrhenius  law.  On  the  other  hand,  for 
reactions  such  as  reactions  5-8  in  Eqs.  (II -3.1)  proceeding  to  the 
right  where  the  rate  constant  does  not  have  an  exponential  term  it 
can  be  expected  that  turbulence  would  effect  the  reaction  rates  in  the 
sense  of  increasing  the  probability  of  useful  collisions.  These 
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statements  impLy  that  all  of  the  reactions  of  interest  in  the  hydrogen- 
oxygen  reaction  except  the  recombination  reactions  (5-8  to  the  right) 
will  be  accelerated  by  turbulence. 

It  would  be  of  fundamental  interest  to  conduct  experiments, 
e.g.  ,  in  nozzles  expanding  the  products  of  hydrogen-oxygen  combustion 
with  various  degrees  of  turbuLence  in  the  stream,  to  determine  the 
effect  of  turbulence  on  reaction  rates  of  practical  importance  in  super¬ 
sonic  combustion. 


SECTION  IV 

EXPERIMENTAL  RESULTS  OF  COMBUSTION 

In  order  to  investigate  experimentally  the  process  of  supersonic 
combustion  controlled  by  mixing  of  hydrogen  and  air,  a  series  of  tests 
is  in  progress  at  the  Aerodynamics  Laboratory  of  the  Polytechnic 
Institute  of  Brooklyn.  Some  of  the  results  of  such  investigation  are 
presented  here. 

EXPERIMENTAL  APPARATUS 

The  tests  presented  here  were  performed  in  an  axially  sym¬ 
metric  jet.  The  hydrogen  was  injected  axially  through  subsonic  jet  at 
the  axis  of  the  air  (see  Figure  IV- 1).  The  air  Mach  number  was  chosen 
equaL  to  1.54  in  order  to  change  sufficiently  some  of  the  parameters 
which  are  important  for  chemical  reaction  and  which  affect  significantly 
the  reaction  time.  The  maximum  air  stagnation  temperature  available 
in  the  facility  was  roughly  1660°K;  therefore,  the  maximum  static 
temperature  of  the  air  couLd  be  changed  for  M=  1.54  between  room 
temperature  and  1110°K.  With  this  Mach  number  it  was  possible  to 
change  the  static  pressure  in  the  mixing  region  from  vaLues  of  the  order 
of  1/20  atm  to  severaL  atms.  However,  for  safety  reasons  the  maximum 
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pressure,  which  permits  optical  observation  through  the  glass  windows, 
was  limited  to  1.3  atm.  The  temperature  of  the  hydrogen  was  also 
variable  between  1110°K  and  324°K. 

The  hydrogen  and  air  jets  discharge  into  a  large  chamber  where 
the  external  pressure  can  be  carefully  controlled.  Therefore,  the 
mixing  and  combustion  process  take  place  in  an  open  jet  region,  which 
presents  several  advantages. 

(a)  It  permits  the  static  pressure  along  the  external 
part  of  the  jet  to  be  kept  constant,  independently 
of  the  combustion  process. 

(b)  It  simplifies  the  optical  observation  of  the  flow. 

The  measurements  made  were  pressure,  velocities  and  temper¬ 
ature  of  the  two  streams.  The  combustion  process  was  observed  by 
direct  photographs  and  schlieren  photographs.  Temperature  survey  was 
also  made.  Direct  photographs  which  are  important  in  understanding 
the  processes  in  the  flow  were  taken  by  a  long  time  exposure  (about  1-7 
sec.)  of  the  flame. 

The  purpose  of  the  tests  was  to  obtain  some  experimental  infor¬ 
mation  on  the  chemical  and  mixing  process.  The  mixing  process  in  the 
tests  is  complex  because  of  the  presence  of  the  boundary  layer  outside 
of  the  hydrogen  injector;  therefore,  initial  tests  were  made  in  order  to 
determine  the  effects  of  the  boundary  layer  on  the  combustion  tests. 

The  presence  of  the  air  boundary  layer  is  important  for  two  reasons; 

(a)  it  increases  the  mixing  process,  and  (b)  the  boundary  layer  air  has 
regions  where  the  static  temperature  is  higher  and  the  velocity  is  lower 
than  the  static  temperature  and  velocity  of  the  air  jet.  It  must  be  noted 
that  in  all  tests  the  walls  of  the  hydrogen  jet  are  preheated  close  to  the 
stagnation  tempe rature  of  the  air. 

The  first  part  of  the  test  program  has  been  performed  at  a  static 
pressure  of  0.9  atm.  For  this  condition  the  chemical  behavior  is 
described  by  the  data  of  Figure  IV- 2,  which  gives  the  induction  and 
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energy  release  times  as  functions  of  local  temperature.  In  the  range 
of  temperatures,  and  for  the  pressure  chosen  for  the  tests,  these  times 
are  practically  independent  of  the  chemical  composition  of  the  mixture 
(value  of  r\);  therefore,  it  is  possible  to  obtain  indications  on  the  process 
of  mixing  and  reaction  also,  if  the  air  and  hydrogen  concentrations  of 
the  mixture  are  not  constant  and  are  not  exactly  known  along  any  given 
streamline  of  the  fLow.  Figure  IV- 2  aLso  shows  that  the  energy  release 
time  is  very  short,  and  the  induction  time  represents  the  larger  portion 
of  the  total  combustion  time.  Rapid  changes  of  reaction  time  take  pLace 
in  the  range  between  800°  and  900°K.  For  example,  a  variation  of 
temperature  between  940°K  and  830°K  changes  the  reaction  time  between 

*.4  _  4 

2.0x10  to  10x10  ,  which  corresponds  to  variations  of  distance  in 

the  flow  of  several  inches.  This  large  variation  of  travel  length  permits 
an  interpretation  of  the  data  obtained  from  photographic  observations. 


Several  series  of  tests  were  performed  for  the  cited  air  Mach 
number  of  1.54  and  for  hydrogen  Mach  number  variable  but  less  than 
one.  In  each  series  of  tests  the  stagnation  temperature  of  the  air  and  of 
the  hydrogen  were  kept  constant,  and  the  velocity  of  the  hydrogen  u.  was 
changed. 


Several  series  of  tests  were  performed  corresponding  to  different 
combinations  of  stagnation  temperature  of  the  air  and  of  the  hydrogen. 

The  mixing  analysis  does  not  permit  accurate  determination  of 
the  mixing  process;  however,  it  does  permit  obtaining  some  relation 
between  velocity  stagnation  temperature  and  static  temperature  as  function 
of  concentration.  If  0g  and  pj  are  the  partial  mass  densities  of  air  and 
hydrogen  at  a  given  position  where  the  mass  density  is  p,  then  the  average 
velocity  and  stagnation  temperature  can  be  expressed  as 


pu=peue*  Pju.  (IV -1) 

where  ug  and  u.  are  the  initial  velocities  of  air  and  hydrogen  stream  and 
the  effects  of  viscous  Losses  are  neglected.  This  approximation  is 
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satisfactory  in  view  of  the  fact  that  the  exchange  of  momentum  between 
the  air  and  hydrogen  stream  represents  only  a  small  percentage  of  the 
momentum  of  the  air  stream,  and  therefore  the  viscous  losses  related 
to  this  exchange,  that  depends  on  the  momentum  exchanged  can  be 
neglected,  and 


H  =  c  T 
o  p  o 


C  c  T  +  p .  c  T 
e  p  s  j  p.  s. 
e  e  }  i 


(IV -2) 


As  it  has  been  discussed  in  Section  III,  these  expressions  are  consistent 

with  the  concept  of  the  mixing  length  introduced  by  Prandtl  when 

molecular  transport  phenomena  can  be  neglected.  Therefore,  if  K  =  - — -  , 

P 


T 

o 


H  +K(H.-H  ) 
e  J  e 


cp  +fc 


Pi 


-c  )K 
P  ' 


(IV -3) 


and 


u  =  u  + 
e 


(u.-u  )K 
J  e' 


(IV -4) 


T  = 


c  T  +K(c  T  -  c  T  )  -  4  us 
Pe  se  Pj  sj  Pe  Se  2 


c  >K  + 

P„  P 


e 


(IV -5) 


F rom  these  relations  and  from  the  relation  between  reaction  time  and 
static  temperature  T  the  effect  of  the  variation  of  H.  and  H  and  of  the 

J  e 

concentration  on  the  flame  front  can  be  investigated. 

In  the  first  series  of  tests  the  stagnation  temperature  of  the 
air  was  chosen  to  be  il70°K  corresponding  to  an  air  velocity  of 
2870  ft/sec  and  to  a  static  temperature  of  834°K.  The  stagnation 
temperature  of  the  hydrogen  was  chosen  equal  to  985°K.  The  Mach 
number  of  the  hydrogen  stream  is  low,  and  the  static  temperature 
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changes  only  slightly  when  u.  changes,  and  is  close  to  the  stagnation 
temperature;  therefore,  in  this  series  of  tests  the  static  temperature 
of  the  hydrogen  is  larger  than  the  static  temperature  of  the  air.  For 
these  conditions  the  reaction  time  was  expected  to  be  of  the  order  of 

-  4 

10  .  During  this  time  the  air  in  the  boundary  layer  has  time  to  travel 

several  inches  and,  therefore,  because  of  the  small  boundary  layer 
thickness,  it  could  be  expected  that  the  boundary  layer  wou'd  be  mixed 
with  the  main  flow  before  combustion  starts.  Some  of  the  results  of 
such  series  of  tests  are  presented  in  the  next  figures. 


Figure  IV-3a  presents  a  direct  photograph  of  the  process  of 
-  ui 

combustion  for  low  velocity  ratio,  =  0.3  .  A  drawing  of  the  flame 
front  for  the  same  condition  is  given  in  Figure  IV-3b.  For  this  con¬ 
dition  the  velocity  of  the  hydrogen  jet  is  equal  to  860  ft/sec.  The  flame 
starts  at  the  axis  at  a  distance  of  3.5"  downstream  of  the  exit  of  the 
hydrogen  injector;  it  propagates  rapidly  across  the  central  region  of  the 
mixing  zone  and  then  slows  down  sharply  when  the  flame  reaches  the 
outside  region  of  the  jet.  For  a  short  period  the  flame  width  increases 
very  slowly.  At  a  distance  of  the  order  of  7"  from  the  exit  of  the  ejector 
the  slope  of  the  flame  front  changes,  and  the  velocity  of  propagation  in 
a  direction  normal  to  the  jet  increases  again.  Downstream  of  this  point 
the  flame  width  grows  approximately  at  a  constant  rate,  and  the  flame 
front  is  inclined  at  an  angle  of  approximately  7°-8°  with  respect  to  the 
axis  of  the  jet. 


The  shape  of  the  flame  is  consistent  with  what  could  be  expected 
from  the  results  of  an  anaLysis.  In  Figure  IV -4  the  flow  travel  required 
for  combustion  is  given  as  a  function  of  local  average  concentration  for 
the  conditions  of  this  test.  This  flow  travel  has  been  calculated  by 
means  of  Eqs.  (IV-4)  and  (IV-5)  on  the  basis  of  the  data  of  Figure  IV-2, 
and  does  not  take  into  account  the  effect  of  the  boundary  layer. 

At  the  center  of  the  flame  the  hydrogen  concentration  is  very 
high;  therefore,  the  static  temperature  is  high  and  the  velocity  is  low. 
The  flow  travel  for  combustion  given  by  Figure  IV-4  here  is  of  the 
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FIG.  IV-3a.  Direct  Photograph  of  Air  Hvdrogen  Combustion 

T  =1170°K;  T  = y»buK;  u./u  =0.3 
s  s.  l'e 

e  J 
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Length  versus  Mass  Concentration  for  Combustion  T  =1170  K;  Tg  =  970 


order  of  1.9".  Therefore,  about  1  6"  of  flow  traveL  are  required  before 
some  air  reaches  the  center  of  the  jet.  When  the  concentration  decreases, 
the  flow  travel  required  for  combustion  increases;  however,  the  mixing 
length  decreases,  Hydrogen  concentrations  of  the  order  of  .50  or  Lower 
are  reached  almost  immediately  downstream  of  the  injector  at  the  outer 
edge  of  the  mixing  region;  therefore,  for  those  conditions  the  totaL  flow 
travel  between  the  end  of  the  ejector  and  the  fLame  front  corresponds 
roughly  to  the  fLow  traveL  required  for  combustion.  The  fLow  traveL 
required  for  combustion  for  K  =  0.50  corresponds  to  4.4".  This  explains 
the  shape  of  the  fLame  in  the  central  region  of  the  mixing.  When  the 
concentration  decreases  the  flow  traveL  for  combustion  becomes  very 
Large;  therefore,  the  veLocity  of  propagation  of  the  fLame  across  the  jet 
decreases  very  rapidly. 

The  region  of  the  flame,  where  the  width  of  the  flame  remains 
practicaLLy  constant,  can  be  attributed  to  the  presence  of  the  boundary 
layer  in  the  air  stream,  generated  along  the  external  surface  of  the 
injector.  The  air  boundary  layer  near  the  walL  of  the  injector  has  larger 
static  temperature  and  lower  velocity  than  the  fLow  outside  of  the  boundary 
layer.  As  a  consequence,  in  the  region  where  the  mixture  contains 
boundary  layer  air,  combustion  times  and  flow  velocities  are  smaller; 
therefore,  flow  travels  for  combustion  are  much  smaller.  Then  the 
flame  front  propagates  fast  into  regions  of  low  concentrations.  When 
the  air  velocity  increases  and  static  temperature  decreases,  the  flow 
traveL  for  combustion  for  a  given  vaLue  of  concentration  increases  and, 
therefore,  the  fLame  front  moves  into  regions  of  higher  concentrations 
than  for  the  case  of  a  mixture  having  boundary  Layer  air.  As  a  con¬ 
sequence,  the  transversal  dimension  of  the  flame  increases  very  sLowLy 
because  the  hydrogen  concentration  must  increase  along  the  flame  front. 
When  the  flame  reaches  a  region  where  the  air  has  free  stream  condi¬ 
tions,  the  flame  front  propagates  along  a  Line  of  roughly  constant  con¬ 
centration,  On  this  basis  the  flame  front  permits  information  to  be 
obtained  on  the  mixing  process.  Mixing  anaLysis  indicates  that  the 
mixing  for  these  conditions  propagates  at  angles  of  the  order  of  7°. 
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When  the  .clccity  o£  the  jet  increases,  then  the  fLame  changes 
shape,  as  shown  in  Figure  IV-5.  The  flame  moves  downstream  of  the 
injector  and  then  propagates  in  similar  fashion  as  before.  Again  the 
shape  of  the  flame,  the  distance  from  the  discharge  section  of  the 
injector,  is  consistent  with  the  mixing  process  and  the  calculated 
reaction  rates. 

In  the  photograph  shown  the  velocity  of  the  hydrogen  u;  is  equal 
i  ui  ^ 

to  1430  ft/ sec  and  0.5 .  The  stagnation  temperature  of  the  air  is 

‘  ■ 

1I90°K,  and  of  the  hydrogen  is  970°K.  The  relation  between  flow 
travel  for  combustion  and  concentration  for  these  conditions  is  given 
also  in  Figure  IV-4.  The  flame  starts  at  the  axis  of  the  mixing  zone 
at  a  distance  of  roughly  7,5".  There  the  travel  for  combustion  corre¬ 
sponds  approximately  to  4.5";  therefore,  the  mixing  travel  is  equal  to 
3".  If  Uj  continues  to  increase,  the  flame  moves  gradually  downstream 
of  the  exit  of  the  injector.  For  high  ■j-j-i  the  mixing  process  becomes 
slower  and  the  combustion  starts  at  the  outer  edge  of  the  mixing  region. 

In  Figure  IV-6  the  distance  x  in  inches  measured  from  the 

u  i 

discharge  section  of  the  injector  is  plotted  as  a  function  of  — d  .  For  the 

e 

case  of  uj  =  ue  the  velocity  of  the  hydrogen  is  2820  ft/sec,  and  the  static 
temperature  of  the  hydrogen  is  890°K.  For  these  conditions  the 
temperature  of  the  mixing  is  also  of  the  order  of  872°K,  and  the  com- 
bustion  will  occur  after  4.5x10  sec.  Then  the  distance  would  be  of 
the  order  of  13".  This  value  does  not  agree  exactly  with  the  experi¬ 
mental  results  shown  in  Figure  IV-6  which  gives  a  distance  of  the 
order  of  10".  The  difference  is  due  to  the  presence  of  the  boundary 
Layer  that  slightly  increases  localLy  the  static  temperature  and 
decreases  the  velocity  of  the  mixture. 

If  the  velocity  u-  increases  to  a  value  of  2u  .  then  the  static 
J  ^ 

temperature  ci  the  hydrogen  becomes  equal  to  834°K  and  the  reaction 

—  4 

time  becomes  o*  the  order  of  6.7x10  ;  therefore,  the  Length  must 

become  of  the  order  of  24"  The  experimental  data  agree  approxi¬ 
mately  with  the  values  of  the  distances  given  by  the  anaLysis. 
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FIG.  IV-5.  Direct  Photograph  of  Air  Hydrogen  Combustion 

T  =11 90°K;  T  =  970°K;u./u  =0.5 
s  s.  re 

e  J  J 
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In  order  to  understand  the  interplay  between  mixing  and  com¬ 
bustion  phenomena  of  hydrogen  and  air,  it  can  be  useful  to  discuss  two 
other  series  of  tests  having  different  stagnation  conditions.  In  the 
second  series  the  stagnation  temperature  of  the  hydrogen  is  the  same 
as  the  static  temperature  of  the  air,  and  this  temperature  is  somewhat 
higher  than  in  the  first  series;  therefore,  the  reaction  rates  are 
faster.  In  the  third  series  of  tests  the  stagnation  temperature  of  the 
hydrogen  is  much  lower  than  the  static  temperature  of  the  air. 

T  =  1 3 20°K  T  =  954°K  -  1000°K  series  2 

s  s . 

e  J 

T  =  1320°K  T  =  362°K  series  3 

s  s . 

e  J 

For  the  second  series,  because  the  static  temperature  of  the 
hydrogen  and  of  the  air  are  about  the  same.  The  reaction  rate  is 

-  4 

constant  and  of  the  order  of  1.7x10  independent  of  concentration. 
Photographs  for  vaLues  of  uj  =  0.28,  0.92,  1.21  ue  are  shown  in 
Figures  IV-7,  -8,  and  -9,  respectively. 

For  low  vaLues  of  Uj  the  combustion  starts  first  at  the  center. 
Figure  IV-7.  The  distance  is  roughly  2.6"  from  the  exit  of  the  jet. 

In  this  test  the  hydrogen  temperature  is  equal  to  1000°K.  and  the 
reaction  time  at  the  axis  is  l.OxlO"4.  The  velocity  is  860  ft/sec; 
therefore,  the  reaction  length  is  of  the  order  of  1";  and  the  mixing 
length  is  of  the  order  of  1.6",  that  is,  equal  to  the  value  determined  in 
the  preceding  test  for  about  the  same  vaLue  of  uj/ue.  The  fLame  shape 
is  as  shown  schematically  in  Figure  IV- 10.  The  flame  propagates 
rapidly  across  the  jet.  The  flame  front  depends  only  from  the  veLocity 
distribution  because  the  temperature  of  the  mixture  is  approximately 
constant  everywhere.  The  velocity  outside  is  3100  ft/sec,  therefore 
the  flame  reaches  the  outside  of  the  mixing  at  a  distance  of  the  order 
of  about  6",  which  corresponds  to  a  reaction  time  of  1,7x10  ,  This 
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FIG.  IV-7.  Direct  Photograph  of  Air  Hydrogen  Combustion 

T  =  I300°K;  T  =970°K;u./u  =0.28 
s  s.  i  e 

e  J 
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FIG.  IV-8.  Direct  Photograph  of  Air  Hydrogen  Combustion 

T  =  1 300°K;  T  =970°K;u./u  =0.92 
s  s.  )'e 
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FIG.  IV -9.  Direct  Photograph  of  Air  Hydrogen  Combustion 

T  =  1300°K;T  =  960°K;  u./u  =1.21 

s  s.  j  e 

e  J 
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agrees  with  the  measured  vaLue.  The  shape  of  the  fLame  downstream  of 
this  point  gives  the  measure  of  the  velocity  of  turbulent  diffusion  of  the 
hydrogen  normally  to  the  stream. 

If  the  velocity  uj  increases,  then  the  time  required  for  initiation 
of  the  combustion  at  the  axis  increases;  however,  the  outside  conditions 
do  not  change.  This  is  shown  from  photographs  in  Figures  IV-8  and  -9. 

If  uj  =  0.92  ue  ,  the  static  temperature,  and  therefore,  the  reaction  time 
are  about  the  same  as  before;  however,  the  reaction  length  at  the  center 
is  Larger  because  of  the  increase  in  velocity.  The  combustion  starts  at 
a  distance  of  the  order  of  6".  The  flame  outside  has  the  same  shape  as 
before. 

For  the  condition  u;=1.21  ue  the  hydrogen  temperature  is 
sLightLy  Lower.  The  flame  has  the  same  general  configuration  as  for 
uj  =  0.92  ue  ;  however,  it  starts  slightly  downstream  of  the  preceding 
case.  The  difference  is  of  the  order  of  l".  This  would  md’ •'■»*■«»  a  varia¬ 
tion  or  mixing  length  of  the  order  of  13%;  this  is  in  qualitative  agree¬ 
ment  with  the  mixing  tests,  ALL  the  photographs  indicate  that  the  veLocity 
of  propagation  of  the  flame  normal  to  the  stream  increases  moving 
downstream  with  the  flow.  This  can  be  attributed  to  a  graduaL  axiaL 
pressure  rise  connected  with  the  heat  reLease  due  to  combustion.  These 
tests  indicate  that  the  reaction  time  calculated  is  in  agreement  with 
experiments.  It  aLso  confirms  the  resuLts  of  the  mixing  experiments 
which  indicate  that  the  mixing  is  fast  when  uj  <  ue  ,  and  decreases  only 
gradually  when  uj  increases  without  any  signiicant  changes  at  uj  =  ue  . 

In  aLL  the  tests  presented  before  the  process  of  heat  reLease 
is  relatively  slow,  and  therefore  the  locaL  pressure  rise  due  to  com¬ 
bustion  is  smaLl  and  graduaL.  The  schlieren  pictures  for  such  fLows 
indicate  that  the  compression  generated  by  combustion  waves  does  not 
form  an  envelope,  and  shocks  are  not  found  in  the  flow.  However,  the 
absence  of  shock  is  directLy  dependent  on  the  mechanism  of  heat  reLease 
and  on  the  dimensions  and  type  of  boundary  of  the  air  jet.  Shocks  can 
be  produced  either  by  changing  the  geometry  of  the  air  jet  or  by 
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accelerating  and  increasing  the  mechanism  of  heat  release.  Photographs 
in  Figures  IV - 1  1  and  -12  are  proof  of  such  possibilities.  Both  photo¬ 
graphs  are  a  composition  obtained  by  superimposing  a  schlieren  photo¬ 
graph  to  a  flame  photograph  taken  from  the  same  experiment. 

The  windows  for  the  two  photographs  are  not  of  the  same  dimen¬ 
sions  and  the  schlieren  axis  is  perpendicular  to  the  axis  of  the  flow, 
while  the  flame  photograph  is  inclined  with  respect  to  the  axis  of  the  flow 
as  shown  in  Figure  IV - 1 .  Therefore  a  difference  of  angle  exists  for  the 
two  data.  However  the  superposition  is  qualitatively  correct  and  gives  a 
correct  impression  of  the  relation  between  heat  release  and  shock 
formation. 

Two  different  tests  have  been  performed  in  order  to  indicate  the 
effect  of  the  two  parameters,  mechanism  of  heat  release  and  boundary 
conditions , 

The  photograph  of  Figure  IV-11  corresponds  to  stagnation  con¬ 
ditions  Ts  =  1440°K,  Tsi=935°K,  and  u:  =  0.28  ue  .  For  these  conditions 
the  reaction  rate  is  about  twice  as  fast  as  for  the  case  of  Figure  IV -7 
and  therefore  the  heat  release  is  faster.  At  the  same  time,  the  mixing 
rate  increases  when  combustion  takes  place  because  of  the  pressure 
rise  due  to  heat  release;  therefore,  as  soon  as  the  combustion  rate 
increases  the  mixing  rate  also  increases,  and  consequently  the  heat 
release  becomes  much  faster.  The  faster  rate  of  heat  release  produces 
stronger  pressure  rise  and,  therefore,  produces  shocks.  For  the  case 
shown  in  Figure  IV-12  the  stagnation  temperature  of  the  air  is  equal  to 
1170°K  and  the  hydrogen  temperature  to  362°KS  therefore,  the  com¬ 
bustion  rate  for  a  given  concentration  is  lower.  However  in  this  test 
the  mixing  process  has  been  strongly  augmented  by  creating  a  positive 
pressure  gradient  along  the  axis  of  the  mixing  stream. 

The  hydrogen  jet  at  the  exit  of  the  injector  is  subsonic,  while 
the  air  jet  is  supersonic,  therefore,  for  a  short  length  the  central- 
part  of  the  mixing  region  is  subsonic.  If  a  pressure  rise  is  produced 
in  the  flow,  in  some  region  not  too  far  from  the  exit  of  the  injector. 
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FIG.  rv-ll. 


Direct  and  SchLieren  Photographs  of  Air  Hydrogen  Combustion 
T  =  1440°K;  T  =  935°K;  u./u  =0.28 

Se  si  y  e 
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FIG.  IV-12.  Direct  and  SchLicren  Photographs  of  Air  Hydrogen  Combustion 

T  =  1280°K;T  =  380°K;u./u  =  0 .27;  p  /p,  =  1 A 

s  s .  re  a  b 
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this  disturbance  can  travel  upstream  through  the  subsonic  central,  part 
and  from  the  center  propagates  into  the  supersonic  region  Then  the 
velocity  of  the  central  part  of  the  mixing  tends  to  decrease  because  of 
the  pressure  rise,  while  the  mixing  tends  to  produce  the  opposite 
effect.  As  a  consequence  the  mixing  is  strongly  increased,  specifi¬ 
cally  if  the  velocity  at  the  center  is  low  I:;  ;b<  present  test  the  static 
pressure  rise  has  been  produced  by  increasing  the  static  pressure  in 
the  discharge  chamber  to  a  value  of  1  A  the  static  pressure  of  the  air 
stream  at  the  exit  of  the  nozzle.  Then  a  shock  is  produced  at  the  end 
of  the  nozzLe  that  travels  toward  the  axis  and  reaches  the  axis  some¬ 
what  downstream  of  the  region  where  combustion  starts.  The  shock 
produces  a  pressure  rise  in  the  central  part  of  the  mixing  where  the 
stream  is  still  subsonic.  The  pressure  rise  decreases  the  velocity  and 
the  mixing  sharply  increases  the  heat  release  While  shocks  would  not 
be  formed  in  this  region  in  absence  of  combustion  the  rapid  com¬ 
bustion  process  produces  sufficient  changes  ir.  density  and  pressure  to 
create  shocks  as  shown  by  the  photograph 

Consider  now  test  series  3.  For  this  scries  the  hydrogen  is 

cold,  therefore,  combustion  can  occur  only  when  the  concentration  is 

low.  The  flow  travel  for  combustion  as  function  of  concentration 

increases  very  rapidly  for  values  of  concentrations  between  0  and  0.1, 

because  this  change  corresponds  to  a  temperature  change  between 

942°K  and  600°K  for  u,  =  0.1  u  and  between  042°K  and  560°K  for 

J  “ 

uj  =  0.78  ue  ,  Therefore  combustion  can  take  place  only  ;-i  regions 
where  the  concentration  is  somewhat  less  than  0  1C  The  velocity  of 
the  jet  for  K  =  0.1  changes  from  2bb7  ft/sec  to  2983  ft/sec  when  uj 
changes  from  0,1  to  0.78  ua  ,  This  change  c!  c.  ;s  small  and  therefore 
it  can  be  assumed  that  the  jet  velocity  ai  K<0.!  is  about  the  same 
independently  of  the  value  of  u-  Figures  IV  ••  i  2  -14.  and  -15  show 

the  flame  shape  for  '  =  0.78.  0  525  and  0  1  respectively.  For  high 
velocity  of  the  hydrogen  jet  (uj/ue  =  0  arlci  0  3231  the  concentration 
below  0.1  is  reached  first  in  the  outside  region  o:  the  mixing  zone,  and 
the  flame  initiates  there.  For  the  temperatures  and  velocities  of  the 
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FIG.  IV- 14.  Direct  Photograph  of  Air  Hydrogen  Combustion 

T  =  1320°K;T  =  340°K;u./u  =0.525 

s  s.  J  e 

e  j 
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FIG.  rv-15. 


Direct  Photograph  of  Air  Hydrogen  Combustion 

T  =  I3  20°K;  T  =380°K;u./u  =0.1 
s  s.  J  e 
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experiments  the  fLow  travel  for  combustion  for  K<0.1,  determined  on 
the  basis  of  the  calculated  reaction  rates,  is  equal  to  6.5",  The  photo¬ 
graphs  of  Figures  IV-13  and  -14  give  for  this  distance  a  value  of  5.5". 
The  smaller  value  of  flow  travel  for  combustion  given  by  the  experi¬ 
ments  is  justified  by  the  presence  of  the  boundary  layer. 


For  very  low  values  of  u-  ,  — *-  =  0.1,  the  combustion  starts  at 

J  ue 

th.e  center.  This  indicates  that  a  very  short  flow  travel  is  required  to 
completely  mix  hydrogen  and  air  when  the  velocity  of  the  hydrogen  is 
small,  so  that  very  low  hydrogen  concentration  is  obtained  at  the  axis. 
Again  this  rapid  mixing  can  be  attributed  to  the  presence  of  the  boundary 
layer  in  the  air  stream.  For  very  low  uj  the  mass  of  the  hydrogen  jet 
is  of  the  same  order  of  the  mass  of  the  air  boundary  layer,  then  the 
mixing  is  fast  and  the  travel  is  short,  because  of  the  lower  velocity  of 
the  boundary  layer  air. 


The  last  experimental  investigation  on  combustion  presented 
here  is  related  to  the  study  of  effects  of  pressure  variation  on  reaction 
rates.  For  this  series  of  tests  the  values  of  uj/ue  and  the  stagnation 
temperatures  of  hydrogen  and  air  have  been  kept  constant  while  the 
static  pressure  has  been  changed.  The  mixing  process  is  only  slightly 
effected  by  the  pressure  changes  because  the  flow  is  turbulent,  there¬ 
fore,  the  variation  of  flame  position  can  be  attributed  to  the  effect  of 
the  variation  of  the  static  pressure  on  the  reaction  rate.  In  Section  II 
an  approximate  equation  has  been  derived  for  the  ignition  delay  as 
function  of  temperature  [Eq.  (11-5.1)].  This  equation  relates  the 
induction  time  to  the  values  of  the  static  pressure  and  temperature. 
The  relation  is  a  function  of  the  values  of  the  chemical  reaction  rates. 
The  experiments  permit  measuring  the  distance  between  the  exit  of  the 
jet  and  the  initiation  of  the  flame;  therefore  it  permits  checking  the 
approximation  of  such  an  equation  and  of  the  reaction  rates  used.  In 
order  to  have  an  accurate  indication  of  the  mixing  length,  the  tests 
were  performed  with  air  at  1440°K  and  hydrogen  cold,  380°K,  so  that 
the  combustion  will  start  only  for  very  low  values  of  concentration 
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F(G.  I V- 1 6 .  Comparison 


(below  K  =  0.05) .  F or  such  conditions  the  velocity  of  the  mixture  is 
equal  to  the  velocity  of  the  air,  and  the  travel  required  for  mixing  is 
negligible,  because  of  the  rapid  mixing  process  and  of  the  low  values 
of  concentration  required  for  combustion. 

An  analysis  of  the  photographs  indicates  that  the  flame 
front  has  a  discontinuity.  This  discontinuity  is  due  to  the  presence  of 
the  boundary  layer  in  the  air  stream.  As  it  has  been  explained  before, 
the  small  amount  of  combustion  that  occurs  before  the  discontinuity  is 
due  to  the  low  velocity  part  of  the  boundary  layer.  The  distance  used 
for  the  determination  of  the  reaction  time  is  the  distance  between  the 
end  section  of  the  hydrogen  injector  and  the  discontinuity  on  the  flame 
front  which  corresponds  to  a  mixture  having  a  velocity  equaL  to  the 
free  stream  veLocity  of  the  air.  Figure  IV-16  shows  a  comparison 
between  the  analytical  results  and  the  combustion  times  determined 
from  the  distances  given  by  the  experiments.  The  comparison  indicates 
that  the  analysis  is  in  good  agreement  with  the  experimental  results. 


SECTION  V 
CONCLUSIONS 

The  theoretical  and  experimental  investigation  reported  herein 
indicates  that  for  the  range  of  conditions  investigated  the  analysis  of 
chemical  reactions  between  hydrogen  and  air  is  sufficiently  accurate  to 
permit  fluid  dynamic  analysis  of  the  flow  field  and  interpretation  of 
experimental  data  connected  with  supersonic  diffusion  flames.  The  fLuid 
dynamic  problem  is  complex  and  requires  a  better  quantitative  under¬ 
standing  of  mixing  problems. 

The  experimental  results  indicated  that  turbulent  mixing  of 
hydrogen  air  stream  is  fast  and  at  present  can  be  predicted,  at  least 
qualitatively,  by  a.nalyses;  however,  variations  to  the  incompressible 
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type  of  expression  of  eddy  viscosity  must  be  introduced.  Such  modifi¬ 
cations  seem  called  on  for  physical  grou"^  The  presence  of  air  boundary 
layer  strongly  affects  ihe  mixing  and  nr  .  >e  taken  into  account. 

On  the  basis  of  these  results  it  appears  that  additional  work  along 
the  lines  indicated  in  the  report  is  required  to  improve  the  analysis  of 
the  fluid  dynamic  aspects  of  supersonic  diffusion  flames;  the  chemical 
aspects  of  the  problem  are  well  understood  for  the  range  of  parameters 
investigated  in  the  experiments  reported  here  and  in  related  work  else¬ 
where.  The  question  of  the  effect  of  turbulence  on  the  chemical  kinetics 
is  still  open. 
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